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ABSTRACT 


This  thesis  provides  the  student  of  combat  modelling 
with  a  computer  program  for  a  relatively  simple  combat 
model  that  can  be  used  in  a  classroom  environment  for  study 
and  analysis.  The  model  is  an  aggregated,  force-on- force 
ground-combat  model  that  uses  Lanchester's  aimed- fire 
equations  for  casualty  assessment.  The  original  version 
of  the  model  was  developed  in  1979  in  a  previous  thesis, 
and  shortfalls  in  the  original  model  have  been  overcome, 
along  with  the  addition  of  several  enrichments  to  provide 
added  user  flexibility.  A  user's  manual  is  provided  to 
facilitate  user  access  to  the  model  from  a  permanent  disk 
in  the  W.  R.  Church  Computer  Center. 


TABLE  OF  CONTENTS 


I.  INTRODUCTION  -  6 

II.  THE  ORIGINAL  MODEL -  9 

A.  GENERAL  DESCRIPTION -  9 

B.  ATTRITION  PROCESS  -  9 

C.  UNIT  LOCATION  AND  MOVEMENT - 11 

D.  BATTLE  TERMINATION - 12 

E.  DETECTION  AND  FIRE  ALLOCATION - 13 

III.  MODEL  CHANGES - 3  5 

A.  GENERAL - 15 

B.  CONDITIONALLY-DETERMINISTIC  ATTRITION  -  15 

C.  UNIT  LOCATION  AND  MOVEMENT - 20 

1.  Straight  Line  Routes - - - — — — - 21 

2.  User  Determined  Routes  ———————  21 

D.  BATTLE  TERMINATION - 23 

E.  ALTERNATE  DEFENSIVE  POSITIONS  -  25 

IV.  FUTURE  ENHANCEMENT - 27 

A.  HETEROGENEOUS  FORCES  -  2? 

B.  AMMUNITION  EXPENDITURE  -  28 

C.  ARTILLERY - 29 

APPENDIX  A:  USER'S  MANUAL  -  31 

APPENDIX  B:  PROGRAM  LISTING  -  49 

APPENDIX  C:  DEFINITION  OF  VARIABLES  IN  COMPUTER  PROGRAM-  77 

BIBLIOGRAPHY - 81 

INITIAL  DISTRIBUTION  LIST - 82 


As  the  art  of  combat  modelling  becomes  more  advanced, 
combat  modelers  are  continually  building  more  and  more 
complicated  models.  To  the  beginning  modeler,  the  ability 
to  understand  how  these  models  operate  is  difficult,  if 
not  impossible.  As  a  student  of  combat  modelling,  I  have 
sensed  a  need  for  a  relatively  simple  model  that  could  be 
easily  studied  in  a  classroom  environment.  It  was  with 
this  in  mind  thst  this  project  was  initiated. 

At  present,  there  seems  to  be  no  simple  combat  model 
available  that  demonstrates  the  basics  of  model  building 
to  the  beginning  student.  In  1979 »  Smoler  (a  student 
at  the  Naval  Postgraduate  School)  attempted  to  build  such 
a  model  for  his  thesis  research  (Smoler ,1979).  His  model 
was  a  deterministic,  force-on-force  computer  model  that 
used  Lanchester's  aimed- fire  equations  for  casualty  ass¬ 
essment.  The  general  scheme  of  his  model  is  shown  in 
Figure  1. 

Although  Smoler’ s  model  was  a  bold  attempt  at  a  simp¬ 
listic  combat  model,  it  did  have  some  problem  areas  that 
warranted  investigation.  Some  of  these  problems  were 
first  discovered  during  a  class  project  in  a  combat  models 
class  (OA  4655)  at  the  Naval  Postgraduate  School.  It  is 
the  purpose  of  this  thesis  to  study  Smoler' s  model  and  to 
undertake  an  enrichment  program  that  will  make  the  model 


FIGURE  1.  GENERAL  SCHEME 


more  realistic  and  flexible,  while  maintaining  it's  trans¬ 
parency  and  simplicity. 

Once  enhanced,  the  model  would  then  be  made  available 
for  use  as  an  instructional  tool  for  combat  modelling 
classes.  This  will  be  accomplished  through  the  develop¬ 
ment  of  a  user's  manual  to  facilitate  the  use  of  the  model 
by  students  by  placing  the  model  on  a  permanent  disk  in 
the  W,  R.  Church  Computer  Center  where  it  will  be  easily 
accessed  by  any  desired  user. 

The  remainder  of  this  paper  will  discuss,  in  detail, 
the  problems  found  with  Smoler's  model  and  the  methods 
used  to  solve  these  problems.  It  will  also  outline  some 
new  features  that  have  been  incorporated  into  the  model 
to  allow  user  flexibility.  Finally,  the  User's  Manual 
(Appendix  A)  will  provide  all  the  required  information 
to  enable  even  the  novice  user  to  utilize  the  model. 


II.  THE  ORIGINAL  MODEL 


A.  GENERAL  DESCRIPTION 

Smoler's  original  model  of  land  combat  is  a  deter¬ 
ministic  model  that  plays  combat  between  two  homogeneous 
forces,  a  blue  force  and  a  red  force.  The  blue  force  is 
comprised  of  three  subunits  in  a  static  defense,  with  each 
subunit  armed  with  three  TOW  antitank  missile  systems. 

The  red  force  is  composed  of  three  subunits  of  three  tanks 
each,  attacking  on  pre-planned  routes.  The  battle  takes 
place  on  the  10  x  10  Km  Fulda  Box  that  has  been  developed 
and  used  in  the  STAR  simulation  model  (Wallace  ,1978).  Since 
the  major  thrust  of  this  paper  is  to  alter  the  original 
model  to  a  form  that  will  be  easily  used  for  classroom 
instruction,  a  brief  discussion  of  the  major  components 
of  the  original  model  will  be  presented,  including  problem 
areas  that  have  motivated  model  changes. 

B.  ATTRITION  PROCESS 

The  attrition  process  in  the  original  model  utilizes 
Lanchester  "aimed-fire"  equations  with  variable  attrition 
coefficients.  The  Lanchester  equations  used  are  fairly 
simple  and  will  be  discussed  later.  However,  the  calcu¬ 
lation  of  the  attrition  coefficients  is  of  more  immediate 
interest.  Smoler  used  the  Bonder-Farrell  model  to  calcu¬ 
late  the  coefficients,  A^,  the  rate  at  which  one  firer 
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of  unit  i  kills  unit  j  targets.  These  A^j's  are  computed 
according  to: 

Aij  •  1/E(V 

where  E(T.  .)  is  the  expected  time  for  one  firer  of  unit  i 

■*“  J 

to  kill  one  target  of  unit  j.  The  E(T^j)  is  computed 
using  the  Bonder-Farrell  formula: 

E(Ti j)  =  ta  +  ti  -  th  +  (th  +  tf)/P(KIH)  + 

(tm  +  tf)/P(h|m)  x  ( ( 1 — P( h I h) )/P(KlH)  + 
P(h!h)  -  P) 

where 

t  =  time  to  acquire  a  target 

a 

t^  =  time  to  fire  first  round  following  acquisition 

th  =  time  to  fire  following  a  hit 

t  =  time  to  fire  following  a  miss 

tj  =  time  of  flight  of  a  round 

P  =  probability  of  a  first  round  hit 
P(hlh)  =  probability  of  a  hit  following  a  hit 
P(hlm)  =  probability  of  a  hit  following  a  miss 
P(KlH)  =  probability  of  a  kill  given  a  hit 
This  formula  holds  for  the  conditions  that  the  hit 
probability  of  any  round  depends  only  on  the  result  of  the 
previous  round  and  no  accumulated  damage  is  considered. 
Smoler  also  assumed  that  for  the  TOW  weapon  system, 

P(KIH)  =  1.0  and  P(hlra)  =  P(h|h)  =  P,  thus  reducing  the 
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equation  to: 
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E(Tij)  =  ta  +  *1  +  tf  +  (tm  +  tf)(l-P)/P 
Smoler  also  assumed  that  for  tanks,  P ( K | H )  =  1.0  and 
tj  =  0,  reducing  the  formula  to: 

E(Tf  .)  =  ta  +  t1  +  tm(l-P)/P(hlm) 

Utilizing  these  equations  to  calculate  the  A^'s,  the 
attrition  during  each  time  step  was  computed  using  the 
Euler-Cauchy  differencing  equations  to  approximate  Lan- 
chester's  fo^c e-on-force  attrition  differential  equations. 

0.  UNIT  LOCATION  AND  MOVEMENT 

Red  and  blue  unit  locations  v;ere  handled  in  two  dif¬ 
ferent  ways  in  the  original  model.  The  blue  locations 
were  left  as  user  inputs,  while  the  red  locations  had 
been  pre-determined  by  the  model  builder  and  could  not  be 
altered  by  the  user.  This  allowed  flexibility  of  defensive 
positions,  but  because  of  the  method  of  determining  movement 
routes  for  the  attackers,  this  flexibility  was  limited. 

The  method  utilized  for  route  determination  was,  at 
best,  unrealistic.  For  each  original  red  location,  a 
straight  west  to  east  route  was  calculated.  Each  route 
was  divided  into  40  meter  intervals,  since  each  red  unit 
was  assumed  to  move  that  distance  during  each  10  second 
time  interval.  This  approximated  an  average  rate  of  move¬ 
ment  of  9  m.p.h.  This  method  of  route  selection,  although 
easy  to  implement,  has  several  significant  shortcomings. 

It  does  not  allow  the  attacking  units  to  utilize  terrain 
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features  during  movement.  Also,  it  permits  attackers  to 
move  over  terrain  that  is,  in  real  life,  impassible. 

D.  BATTLE  TERMINATION 

In  any  combat  model,  adequate  battle  termination  rules 
must  be  considered.  Smoler  utilized  two  criteria  to 
terminate  the  battle.  The  first  of  these  was  annihilation 
of  one  of  the  two  forces.  This  criteria  is  reasonable  in 
an  expected  value  model  like  this  and  was  adequately 
handled.  This  is  not  the  case  for  the  second  termination 
rule. 

The  second  rule  for  termination  is  that  the  distance 
between  red  and  blue  forces  becomes  too  small.  The  use 
of  this  criteria  for  a  stopping  rule  is  fine,  provided  it 
is  implemented  properly.  In  the  original  model,  a  geographic 
center  was  calculated  for  each  force.  It  was  these  center 
of  mass  points  that  were  compared  to  determine  if  the 
units  were  too  close.  This  method  can  lead  to  problems, 
such  as  those  that  were  encountered  during  some  initial 
trial  runs.  In  one  battle,  the  red  units,  during  their 
advance,  had  eliminated  one  entire  blue  unit  and  had  suf¬ 
fered  the  loss  of  one  of  their  own  units.  Since  the  red 
units  are  allowed  to  advance  as  much  as  150  meters  apart 
and  blue  units  were  set  up  so  they  were  approximately 
1300  meters  apart,  there  was  considerable  distance  between 
units  on  the  battlefield.  Since  the  two  units  had  been 
destroyed,  the  center  of  mass  computation  allowed  the  red 
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units  to  move  completely  through  the  blue  defense,  while 
both  sides  still  had  forces  remaining  to  fight  with.  This 
reflects  a  problem  with  the  center  of  mass  computation  as 
a  method  for  determining  distances  between  units. 

E.  DETECTION  AND  FIFE  ALLOCATION 

Both  detection  and  fire  allocation  processes  are  hand¬ 
led  well  in  the  original  model  and  no  changes  have  been 
made.  A  detailed  description  of  these  methods  is  contained 
in  Smoler's  Thesis (Smoler, 1979) •  However,  a  brief  and 
general  discussion  of  both  processes  will  be  included  here. 

The  detection  phenomena  is  modeled  in  two  ways.  First, 
a  non-firing  detection  can  occur  as  a  result  of  an  observer’s 
random  search  within  his  designated  section  of  responsibility. 
A  30°  field  of  vision  for  an  observer  is  assumed,  and  the 
probability  that  an  observer  is  looking  in  the  direction 
which  enables  him  to  detect  a  target  is  computed  by  inte¬ 
grating  the  Limicon  Function  over  limits  that  are  ±  15° 
from  the  primary  direction  the  observer  is  looking.  The 
Limicon  Function,  f(0),  is  the  following  probability  den¬ 
sity  function: 

f(6)  =  A  +  B  x  cos(0) 

where 

D  =  assigned  sector  width/2 
A  =  -B  x  cos(D) 

B  =  1  /  2(sin(D)-  D  x  cos(D)) 

0  =  primary  direction  observer  is  looking 


Also,  A  and  B  are  chosen  such  that 


/D 

/  f(e)de  =  1 

-D 

The  second  method  of  detection  is  a  firing  detection. 
This  phenomena  occurs  when  the  following  happens.  If  a 
firer's  location  is  within  ±  15°  of  an  observer's  primary 
direction  of  observation  when  he  is  firing,  he  is  assumed 
to  be  detected  and  is  added  to  the  observer's  target  list. 
This  models  the  detection  of  a  firer  by  locating  a  firing 
signature  of  a  weapon. 

The  fire  allocation  process  is  also  modeled  well  by 
Smoler.  Since  each  firing  unit  is  not  restricted  to  firing 
at  one  target,  a  fraction  of  each  firing  unit  is  allocated 
to  each  target  on  the  firer's  target  list.  This  fraction 
is  determined  as  a  function  of  range  to  each  target  and  of 
predetermined  priorities  of  fire  that  are  set  so  each  unit's 
fire  priority  goes  to  those  targets  to  his  immediate  front. 
Targets  to  the  firer's  flanks  are  then  allocated  a  smaller 
percent  of  the  available  firepower.  Again,  the  detection 
and  fire  allocation  processes  have  been  handled  well  and 
no  changes  to  these  routines  will  be  made. 
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III.  MODEL  CHANGES 


A.  GENERAL 

As  the  previous  section  outlines,  there  were  several 
major  problems  discovered  in  Smoler's  original  model.  In 
order  to  make  the  model  more  flexible  and  classroom  useable, 
several  major  changes  were  found  to  be  necessary.  This 
chapter  will  provide  a  detailed  description  of  these 
changes,  which  include  conditionally-deterministic  attrition, 
unit  locations  and  movement,  battle  termination  and  move¬ 
ment  of  defenders. 

B.  CONDITIONALLY-DETERMINISTIC  ATTRITION 

In  order  to  introduce  a  stochastic  process  into  the 
attrition  computation  for  this  model,  several  options  were 
available.  These  include  the  use  of  a  Markov-process  to 
determine  casualties  or  the  use  of  random  attrition-rate 
coefficients.  After  considering  both  options,  and  in 
keeping  with  the  transparent  nature  of  this  model,  it  was 
decided  that  stochastic  attrition  coefficients  should  be 
utilized,  because  the  attrition-rate  coefficient  is  a 
random  quantity  measuring  a  unit's  fighting  ability  and 
can  be  realized  before  any  given  battle.  The  following 
procedure  was  used  to  develop  a  method  for  random  attrition- 
rate  coefficients. 

The  attrition-rate  coefficient,  A^,  is  a  measure  of 
the  rate  a  firer  of  unit  i  attrits  a  target  of  unit  j. 
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This  also  can  be  interpreted  as  a  measure  of  the  fighting 
ability  of  an  element  of  type  i.  It  is  intuitively  obvious 
that  this  is  a  variable  quantity  that  can  be  affected  by 
many  different  factors,  including  weather,  espirit-de-corps, 
previous  engagements  and  leadership,  to  mention  only  a  few. 
In  a  homogeneous  model,  one  can  attempt  to  capture  these 
random  effects  by  developing  a  distribution  of  initial 
fighting  unit  capabilities,  i.e.,  initial  A.  . 's  for  each 

^  J 

unit.  Since  no  data  is  readily  available  that  captures  this 
phenomena,  an  attempt  was  made  to  fit  a  distribution  to 
an  intuitive  feeling  as  to  how  this  fighting  ability  varies 
from  one  unit  to  the  next  prior  to  a  given  engagement.  It 
was  the  author's  intuitive  feeling  that  the  A-.'s  should  be 
distributed  between  .3  and  .8  with  the  majority  of  the  units 
being  rated  between  .3  and  .6. 

Utilizing  the  above  intuition,  an  attempt  was  made  to 
"fit"  a  distribution  to  these  assumed  fighting  levels. 
Graphically,  the  distribution  would  look  like  Figure  2. 


Figure  2.  Distribution  of  Initial  A^'s 


Initially,  an  attempt  was  made  to  fit  a  Beta  distribution 
to  this  curve.  The  fit  was  fairly  accurate  in  the  middle 
range,  but  was  not  satisfactory  in  the  tails.  The  same 
was  true  for  a  Normal  approximation.  It  was  finally 
determined  that  the  best  fit  was  a  straight  quadratic  fit, 
where  a  Uniform(0,l)  input  variable  could  be  used  and  the 
output  would  yield  the  desired  random  number.  The  quad¬ 
ratic  fit  that  was  used  is: 

Ajj  =  -2U2  +  2U  +  .3 

Once  a  distribution  for  the  attrition-rate  coefficients 
had  been  derived,  the  implementation  into  the  model  was 
accomplished.  Since  it  was  assumed  that  the  fighting 
ability  of  each  unit  is  a  random  quantity  prior  to  a  given 
battle,  it  was  only  necessary  to  obtain  a  realization  of 
the  random  variable  for  each  unit  prior  to  the  initiation 
of  the  battle.  This  realization,  A°^,  is  determined  by 
using  a  random  Uniform(0,l)  number  and  the  above  formula. 

A  new  attrition-rate  coefficient  for  each  unit  is  then 
computed  during  each  time  step  using  the  equation: 

f aJj  (l-r/re)2  for  Oir<re 

Au  =  I 

l  0  for  re«r 

where 

r  =  maximum  effective  range  of  a  firer's  weapon 
r  =  current  range  between  firer  and  target 
A^  =  realization  of  unit's  fighting  ability 


This  equation  is  utilized  because  it  varies  with  range, 
but  it  also  is  a  function  of  . ,  thus  creating  a  different 
attrition-rate  curve  for  each  unit,  depending  on  that  unit’s 
fighting  ability  prior  to  the  battle.  Graphically,  as  an 
example,  the  attrition-rate  coefficient  curve  for  an 
of  .6  and  a  maximum  effective  range,  r  ,  of  3000  meters 
would  look  as  in  Figure  3« 


range (meters) 

Figure  3«  Attrition-rate  coefficient  curve  for 
A?^  =  .6  and  rg  =  3000  m. 

Once  the  above  method  for  determining  attrition-rate 
coefficients  had  been  selected,  it  was  coded  and  included 
as  a  user  option  in  the  program.  Once  the  code  had  been 
implemented,  several  runs  were  made  utilizing  different 
random  number  seeds  to  compare  battle  results  using  this 
method  for  attrition-rate  coefficient  determination.  These 
runs  were  then  compared  with  a  run  using  the  original 
deterministic  method.  The  results  are  summarized  below. 

First,  the  model  was  run  using  the  deterministic  attrit¬ 
ion  method.  For  purposes  of  comparison,  all  locations, 
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movement  rates,  force  levels  and  other  input  variables  were 
held  constant  for  all  runs.  Combatant  attrition  using  the 
deterministic  process  as  a  function  of  time  is  shown  in 
Figure  4. 

Force 
Levels 

Time 

Figure  4,  Attrition  using  deterministic  attrition 

In  this  battle,  termination  occurred  as  a  result  of  forces 
being  too  close  together,  with  force  levels  at  termination 
reflected  in  Figure  4.  This  was  then  compared  to  the  runs 
that  utilized  the  new  attrition  module.  The  results  of 
four  of  these  runs  are  in  Figure  5.  The  A°j's  are  the 
realizations  of  the  random  variable  denoting  a  unit's 
initial  fighting  capability  prior  to  the  battle.  The  three 
numbers  under  the  red  and  blue  headings  are  the  realizations 
for  each  subunit  in  the  battle  on  which  casuality  assessments 
are  computed. 

These  results  show,  in  fact,  that  a  randomly  selected 
does  have  an  effect  on  the  final  outcome  of  the  battle. 

This  module  now  provides  the  user  the  option  of  selecting 
a  eonditionally-deterministic  process  that  still  uses 
Lanchester's  equations  as  the  basis  for  the  attrition 
computation. 
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.33  .77  .46 

m 
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.39  .67  .58 

.63  .78  .36 

KB 

0 

330 

Blue  wins 

.78  .67  .80 

.74  .68  .72 

3.2 

0 

320 

Blue  wins 

.32  .43  .68 

.51  .77  .65 

1.2 

3.8 

640 

Too  close 

Deterministic 

IB 

2.8 

460 

Too  close 

Figure  5.  Comparison  of  Stochastic  and  Deterministic 
Battles 

C.  UNIT  LOCATION  AND  MOVEMENT 

As  was  pointed  out  in  an  earlier  section,  the  unit 
locations  and  movement  logic  that  Smoler  used  tended  to 
be  rigid  and  inflexible.  The  changes  that  have  been  made 
to  rectify  this  and  allow  more  flexibility  for  the  user 
will  be  discussed  here. 

First,  the  user  is  now  responsible  for  inputing  all 
combatant  locations,  both  attackers  and  defenders.  The 
format  for  these  inputs  will  be  specified  in  Appendix  A, 
the  User's  Manual.  Aditionally,  the  user  is  given  the 
option  of  selecting  routes  for  the  red  advance  or  allowing 
the  model  to  calculate  straight  west  to  east  routes  as  in 
the  past.  To  add  more  flexibility,  the  user  will  also  be 
required  to  input  the  rate  of  advance  (vehicle  speed)  for 
the  attacking  forces.  The  inputed  rate  of  advance  is 
used  in  both  methods  of  route  calculation,  thus  controlling 
the* red  attack. 
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The  option  to  use  straight  line  routes  from  west 
to  east  is  handled  in  the  same  manner  as  before.  The  route 
is  broken  down  into  discrete  distance  intervals  from  the 
initial  red  location  straight  to  the  east.  The  only  dif¬ 
ference  is  that  the  length  of  each  interval  is  determined 
by  the  vehicle  speed  that  the  user  inputs.  The  route  for 
each  red  unit  consists  of  125  equal  length  intervals. 

2.  User  Determined  Routes 

The  option  to  allow  the  user  to  select  routes  for 
the  red  advance  has  been  added  to  enable  the  user  to  make 
use  of  available  terrain  and  to  add  realism  to  the  model. 

The  method  for  calculating  these  routes  is  straightforward 
and  is  discussed  next. 

The  user  is  required  to  input  the  original  location 
of  each  red  subunit  and  from  one  to  ten  nodes  he  wishes 
each  attacking  subunit  to  move  through.  This  information, 
along  with  vehicle  speed,  is  used  to  calculate  route  inter¬ 
vals  that  move  the  attacking  unit  through  each  of  the 
designated  nodes.  A  complete  route  would  look  like  the  one 
depicted  in  Figure  6.  The  method  used  to  compute  the  routes 
is  as  follows. 

The  straight  line  ground  distance  between  the  first 
two  adjacent  nodes,  DIST,  is  calculated  as  shown  in  Figure 
7.  The  angle  between  the  desired  direction  of  movement 
and  straight  west  to  east  movement,  a,  is  then  calculated. 
Utilizing  these  quantities  and  the  distance  desired  to 


User  Determined  Routes 


move  during  each  time  step,  DST,  the  distance  to  be  moved 
in  the  x  and  y  direction,  XLN  and  YU!,  is  now  computed  as 
shown  in  Figure  7«  These  distances  are  then  added  to  the 
coordinates  of  the  previous  interval  endpoint,  point  C 
in  the  figure,  to  determine  the  coordinates  of  the  next 
interval  endpoint,  point  D.  This  sane  distance  is  again 
added  to  compute  the  coordinates  of  the  next  endpoint, 
point  E.  This  process  is  continued  until  the  distance 
from  the  last  endpoint  computed  to  the  next  node  is  less 
than  DST.  This  process  is  repeated  between  the  next  two 
nodes  until  the  entire  route  is  completed.  In  order  to 
insure  all  intervals  are  of  equal  length,  the  computation 
of  the  first  interval  between  any  two  nodes  must  be  con¬ 
sidered  separately,  by  taking  into  account  the  distance 
left  over  from  the  last  computation  between  the  previous 
two  nodes.  To  do  this,  the  first  interval  takes  the  rem¬ 
aining  distance,  e,  and  adds  it  to  an  interval  length  of 
DST-e  for  the  first  interval  between  any  two  nodes.  This 
insures  that  each  interval  along  the  route  is  of  length  DST, 
which  is  the  required  length. 

D.  BATTLE  TERMINATION 

As  previously  stated,  the  battle  is  terminated  either 
by  annihilation  of  one  of  the  two  forces  or  by  forces 
getting  too  close.  The  annihilation  criteria  has  not  been 
changed.  However,  to  insure  the  distance  criteria  is 
effective,  two  changes  have  been  made. 

First,  the  distance  between  units  is  no  longer  calculated 
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NODE  COORDINATES 

A  XIC( N- 1 ) ,YIC(  N-  1  ) 

B  XLOC(I.J) 

C  XIC(N) ,YIC(N) 

D  XIC ( N+ 1 ) , YIC ( N+ 1 ) 

E  XIC(N+2),YIC(N+2) 

F  XLOC(I ,  J+  1  ) 


USE 

(N-1)st  Interval  Endpoint 

User  Inputed  Node 

Nth  Interval  Endpoint 
(II+1)St  Interval  Endpoint 

^I+2)nd  Interval  Endpoint 
User  Inputed  Node 


DIST  =  NXL2  +  YL2 
a  =  Tan-1  (Y/XL),  where  Y  =  Iyl| 
e  =  distance  less  than  DST  at  end  of  calculation 
of  intervals  between  adjacent  nodes. 

YLN  =  DST  x  SIN(a)  XLN  =  DST  x  COS(a) 


Figure  7.  Route  Computation 
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as  a  center  of  mass  distance.  Instead,  a  distance  is 
computed  between  each  attacking  subunit  on  v/hicn  casualties 
are  being  computed  that  is  still  alive  and  each  defending 
subunit  that  is  still  in  the  battle.  If  any  of  these 
distances  becomes  too  small,  the  battle  is  considered  to 
have  moved  to  a  "close-in"  combat  (hand-to-hand)  that  this 
program  does  not  currently  model.  For  this  reason,  the 
battle  is  then  terminated.  However,  to  insure  the  attackers 
do  not  pass  by  the  flanks  of  remaining  defenders  and  remain 
outside  termination  distance,  a  check  of  x  coordinates 
for  each  subunit  is  also  made.  If  any  attacking  subunit's 
x  coordinate  places  him  past  the  location  of  the  forward 
most  blue  defensive  subunit  still  in  the  battle,  the  battle 
is  also  terminated. 

The  criteria  for  being  too  close  has  been  left  as  a 
user  determination.  It  is  one  of  the  inputs  that  is 
described  in  the  User's  Manual,  Appendix  A.  This  allows 
for  flexibility  of  breakpoint  distances  which  also  lends 
itself  to  the  study  of  optimum  breaking  distances  for 
various  weapon  systems  on  the  battlefield. 

E.  ALTERNATE  DEFENSIVE  POSITIONS 

The  last  addition  to  the  model,  added  to  increase  the 
flexibility  and  realism  of  the  battle,  is  the  option  of 
having  the  defending  units  move  to  alternate  positions  if 
the  attackers  close  to  within  a  user  specified  range. 

The  move  is  handled  in  a  simple  and  transparent  manner. 

If  the  user  selects  this  option,  the  defender  will  move 
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when  the  attacking  forces  close  within  the  breakpoint 
distance,  '.’/hen  this  occurs,  each  defensive  unit  that  is 
still  alive  will  move  to  an  alternate  position  the  user 
has  selected.  The  duration  of  the  move  is  also  a  user 
input.  He  simply  specifies  the  number  of  10  second  time 
steps  he  wants  to  allow  the  blue  forces  to  move.  At  the 
completion  of  the  move,  during  which  the  red  forces  continue 
to  attack,  the  battle  will  continue  until  one  of  the  term¬ 
ination  criteria  discussed  above  has  been  met. 


IV.  FUTURE  ENHANCEMENT 

Although  the  above  changes  have  been  implemented,  and 
a  more  useable  and  flexible  model  has  been  created,  there 
are  several  areas  that  could  still  be  considered  in  future 
work.  These  include  the  introduction  of  heterogeneous 
forces,  explicit  computation  of  ammunition  expenditure 
and  artillery  (indirect  fire).  Possible  methods  for 
employing  these  ideas  are  discussed  below. 

A.  HETEROGENEOUS  FORCES 

The  current  model  involves  combat  between  homogeneous 
forces  only.  In  other  words,  each  opposing  force  is 
comprised  entirely  of  one  weapon  system  type.  Added 
flexibility  could  be  attained  by  allowing  heterogeneous 
force  structures  on  both  sides.  This  would  enable  the  user 
to  investigate  the  effect  of  different  force  mixes  on 
battle  outcome. 

The  introduction  of  different  weapon  systems  within 
a  single  unit  would  require  extensive  changes  to  the  att¬ 
rition  process.  Although  Lanchester  equations  for  aimed- 
fire  could  still  be  utilized,  casualty  assesment  against 
each  weapon  system  type  by  each  opposing  weapon  system 
type  would  have  to  be  calculated.  The  total  attrition  of 
any  particular  unit  would  then  be  the  sum  of  the  damage 
assessed  to  each  weapon  system  of  that  unit.  Aditionally, 
the  force  level  for  each  weapon  system  would  be  required 


as  an  input  as  well  as  a  separate  set  of  hit  and  kill 
probabilities  for  each  weapon  system  type  against  each 
type  of  target. 

Introduction  of  this  feature  would  create  a  more  comp¬ 
licated  model,  adding  realism  but  detracting  from  the  current 
simple  and  transparent  form.  Since  the  purpose  of  the 
current  effort  has  been  to  maintain  this  simplicity,  this 
option  has  not  been  included,  but  could  be  considered  for 
future  work. 

B.  AMMUNITION  EXPENDITURE 

The  combat  process  is  a  complicated  and  intricate 
process.  One  of  the  hardest,  and  not  many  times  attempted, 
areas  to  model  is  the  logistic  area.  However,  in  the 
model  currently  being  studied,  the  problem  of  ammunition 
expenditure  could  be  modeled. 

Since  the  model  is  an  expected  value  type  model,  it 
would  not  be  effective  to  model  the  amount  of  ammunition 
expended  by  simply  counting  bullets  as  they  are  fired. 
Actually,  in  this  type  of  model  it  would  not  be  possible 
to  count  each  round,  because  of  the  aggregated  nature  of 
each  unit.  The  only  way  to  model  ammunition  useage  is  to 
model  it  in  the  same  way  as  attrition.  This  simply  means 
that  each  unit  would  have  to  be  given  a  starting  level  of 
ammunition  on  hand,  a  basic  load,  and  each  time  the  unit 
fired,  the  expected  ammount  of  ammunition  expended  would 
be  subtracted  from  what  is  on  hand.  The  ammount  expended 
would  necessarily  be  a  function  of  the  size  of  the  unit 


firing,  the  number  of  targets  fired  at  and  the  rate  this 
unit  is  firing.  This  can  be  directly  correlated  to  the 
attrition-rate  coefficient  for  the  firing  unit,  '.'/hen 
ammunition  on  hand  reaches  the  zero  level,  the  unit  would 
have  to  be  removed  from  the  battle  permanently  or  for  some 
specified  time  period,  to  simulate  resupply. 

C.  ARTILLERY 

In  any  armed  conflict,  there  are  more  ways  to  inflict 
casualties  than  just  direct  fire.  Many  types  of  indirect 
fire  are  utilized  on  the  battlefield,  including  artillery, 
close  air  support  and  naval  gunfire.  To  include  this  type 
of  play  in  this  model,  primarily  artillery,  would  involve 
a  major  effort  with  large  program  additions. 

To  model  artillery  in  an  expected  value  model  that 
already  uses  Lanchester  equations  for  direct  fire  would 
suggest  the  use  of  Lanchester* s  area- fire  equations.  These 
equations  are  no  more  complicated  to  handle  than  the  airaed- 
fire  equations  that  are  already  in  use.  However,  other 
model  considerations  would  have  to  be  investigated.  These 
would  include  location  of  artillery  units,  whether  to  use 
forward  observers  for  target  location  or  only  use  pre¬ 
planned  strikes  or  both  and  fire  allocation  procedures 
to  be  used.  Also,  a  significant  ammount  of  data  would  need 
to  be  collected  concerning  weapon  types  to  be  used,  eff¬ 
ective  ranges,  killing  radii  and  accuracy  data. 


To  model  artillery  would  be  a  significant  addition  to 
the  model.  However,  this  too  would  effect  the  simplicity 
and  transparency  the  model  currently  possesses.  This 
addition  has  also  been  left  for  future  consideration. 


APPENDIX  A 


User's  Manual 


I.  INTRODUCTION 

Smoler's  model  of  land  combat  is  a  force-on-force 
combat  model  that  utilizes  Lanchester's  aimed  fire  equations 
for  casualty  assessment.  The  battle  is  simulated  on  a 
10  x  10  KM  piece  of  terrain  representing  an  area  east  of 
Fulda,  V/est  Germany.  A  portion  of  the  terrain  map  is 
inclosed  as  Figure  1.  It  is  a  computer  model  that  is 
coded  in  FORTRAN  and  is  available  for  use  in  the  W.  R. 

Church  Computer  Center. 

The  purpose  of  this  manual  is  to  familiarize  the  user 
with  the  model  and  to  provide  the  required  Job  Control 
Language  (JCL)  and  user  inputs  to  run  the  model.  A 
sample  listing  of  program  output  will  also  be  provided  to 
familiarize  the  user  with  the  expected  program  output. 


31 


Model  Terrain  Map 


II.  PROGRAM  STRUCTURE 


Smoler's  enhanced  model  is  a  computerized  model,  coded 
in  FORTRAN,  containing  a  main  program  and  nine  subroutines. 
To  assist  the  user  in  understanding  how  the  model  operates, 
a  brief  description  of  the  function  of  each  of  these  major 
parts  of  the  program  will  be  included. 

A.  MAIN  PROGRAM 

The  main  program  has  several  important  functions. 

All  of  the  input  and  output  functions,  except  the  line-of- 
sight  data,  are  contained  in  the  main  program.  Aditionally, 
the  main  program  is  used  to  structure  all  of  the  other 
functions  during  each  10  second  time  step.  Attrition, 
detection,  movement  and  fire  allocation  are  also  handled 
in  the  main  program.  The  nine  subroutines  provide  needed 
input  numbers  for  the  above  calculations.  The  general  flow 
of  the  program  is  depicted  in  the  flow  chart  shown  at 
Figure  2. 

B.  SUBROUTINES 

There  are  nine  subroutines  included  in  the  model.  Each 
of  these  subroutines  performs  a  distinct  function,  each  of 
which  will  be  discussed  below. 

1.  Subroutine  SETUP 

This  subroutine  is  used  to  read  in  the  terrain 
data  for  the  Fulda  Map.  This  terrain  data  will  be  used 
when  computing  line-of-sight  between  targets  and  observers, 
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FIGURE  2.  FLOW  CIIAPT 


as  well  as  providing  a  grid  system  for  unit  locations  and 
movement. 

2.  Subroutine  LOS 

This  module,  developed  by  Professor  Jim  Hartman 
of  the  Naval  Postgraduate  School,  computes  a  percent  of  a 
target  visible  to  a  particular  observer,  given  the  coor¬ 
dinates  (location)  of  both.  This  visible  fraction  is  used 
in  the  detection  and  attrition  modules  in  the  model. 

3.  Subroutine  KCVEP 

This  subroutine  is  used  by  subroutine  LCS  in  deter¬ 
mining  what  portion  of  a  particular  target  is  covered  by 
th'e  terrain  between  the  target  and  observer.  This  number 
is  used  in  the  detection  and  attrition  modules. 

4.  Subroutine  ETK 

This  module  computes  the  expected  time  for  a  given 
firer  to  kill  a  given  target.  The  calculation  is  a  function 
of  range,  time  of  flight  for  a  round  and  hit  and  kill 
probabilities  for  the  firing  weapon  system.  It  is  a  number 
that  is  used  in  the  computation  of  the  deterministic  attrition 
coefficients,  A.,. 

5.  Subroutine  STOCK 

This  is  used  when  a  user  has  selected  a  stochastic 
attrition  option  to  compute  the  attrition  coefficients 
during  each  time  step.  The  calculation  is  a  function  of 
the  original  stochastically  determined  Ai .  as  well  as  a 
function  of  range. 


7.  Subroutine  LAMBDA 

This  subroutine  is  used  in  conjunction  with  the 
LOS  routine  to  compute  the  detection  rate  of  a  target  by 
an  observer  given  the  percent  of  target  visible  to  the 
observer. 

8.  Subroutine  ROUTS 

This  is  used  to  compute  the  route  of  each  attacking 
red  unit  when  the  user  has  selected  the  option  of  inputing 
attacker  routes.  It  calculates  the  coordinates  of  each 
interval  endpoint  along  the  route,  making  each  interval 
length  (distance  moved  during  a  10  second  time  step)  the 
sane.  The  interval  length  is  determined  by  the  speed  the 
user  has  selected  and  inputed  for  the  current  battle. 

9.  Subroutine  SIEV 

A  subroutine  that  is  used  to  calculate  the  terrain 
elevation  for  a  given  set  of  X,  Y  coordinates.  This  functi 
is  used  in  conjunction  with  the  LCS  subroutine  in  computing 
line-of-sight  between  observer  and  target. 


III.  AVAILABLE  OPTIONS 


The  enhanced  version  of  this  model  has  been  written 
to  allow  for  maximum  user  flexibility  while  maintaining 
the  simple  and  transparent  nature  of  the  model.  To  allow 
this  flexibility,  there  have  been  several  user  options 
incorporated  into  the  model.  Each  of  these  options, 
including  user  responsibilities,  will  be  discussed  here, 
with  the  required  input  data  for  each  being  outlined  in 
the  next  section. 

'AT-  STOCHASTIC  VS  DETERMINISTIC  ATTRITION 

The  user  is  required  to  specify  whether  he  wants 
to  use  stochastic  or  deterministic  attrition  calculations. 
Both  methods  utilize  Lanchester’s  aimed- fire  equations, 
the  difference  being  the  method  of  calculating  the  attrition 
rate  coefficients.  The  deterministic  procedure  uses  the 
Eonder-Farrell  method  of  calculating  the  A .  . ,  while  the 
stochastic  method  selects  an  initial  random  A .  -  for  each 
unit,  and  uses  this  as  a  function  of  range  to  calculate 
the  A^  for  each  time  step. 

E.  ATTACK  ROUTES 

The  second  major  option  available  to  the  user  is 
the  method  of  route  computation  for  the  attacking  forces. 

The  user  has  to  decide  whether  to  allow  the  program  to 
compute  straight  west  to  east  routes  or  to  input  the  routes 
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he  desires  each  attacking  unit  to  follow.  To  select  his 
own  routes,  a  user  must  input  the  number  of  nodes  he 
wishes  to  have  on  each  route  and  the  coordinates  of  each 
of  these  nodes.  The  program  will  then  compute  routes 
through  each  node.  The  nodes  must  be  inputed  in  order 
from  west  to  east  and  should  not  create  an  angle  between 
the  west  to  east  axis  and  the  route  direction  that  exceeds 
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C.  ALTERNATE  DEFENSIVE  POSITIONS 

The  third  option  the  user  must  consider  is  the 
use  of  alternative  defensive  locations.  This  option  permits 
the  user,  if  he  desires,  the  capability  to  move  the  defen¬ 
ders  to  alternate  positions  if  the  attackers  close  within 
some  specified  distance.  This  breakpoint  distance  is 
decided  and  inputed  by  the  user  and  is  also  used  as  the 
distance  for  battle  termination.  The  alternative  to  moving 
the  defenders  is  to  terminate  the  battle  when  the  breakpoint 
distance  is  reached. 

D.  OTHER  INPUTS 

There  are  other  inputs  that  are  required  by  the 
program  that  the  user  must  provide.  These  include  force 
sizes,  weapon  characteristics,  unit  locations  and  hit  and 
kill  probabilities.  The  required  formats  for  all  inputs 
are  outlined  in  the  next  section. 
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IV.  INPUT  DATA 


All  user  input  data  is  read  in  from  cards  at  the  start 
of  the  main  program,  MAIN.  A  brief  description  of  each 
data  element  as  well  as  the  required  format  for  inputing 
this  data  is  discussed  next.  A  sample  data  deck  has 
also  been  included  as  Figure  3,  with  the  referenced  data 
element  number  appended  on  the  right.  This  sample  data 
deck  reflects  the  use  of  all  data  elements.  However, 
not  all  the  elements  described  below  are  required,  and 
those  that  are  needed  only  if  a  particular  option  is 
selected  are  noted  with  an  asteric  next  to  the  data  element 
number. 


Data  Element 


Description  and  Format _ 

Selection  of  attrition  option  and  random 
number  seed.  A  1  is  entered  for  deter¬ 


ministic  or  a  0  for  stochastic  attrition. 


This  number  is  followed  by  a  five  digit 
random  number  seed.  Format  for  this  card 
is  II, IX, 15. 

2.  Number  of  blue  units  (1'IEU)  and  number 

of  red  units  (MRU).  Eoth  of  these 
elements  are  two  digit  integers.  Card 
format  is  l2,iX,I2. 

3«  Effective  weapon  ranges  listed  in  the 

following  order.  Red  minimum,  red 
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maximum,  blue  minimum  and  blue  maximum. 
Card  format  is  4(F6.1,1X). 

Red  force  levels  for  each  unit. 

Card  format  is  F3«l>  IX.  The  card  will 
contain  the  number  of  entries  that  equals 
the  number  of  red  units  (XRU)  in  data 
element  2. 

Type  of  route  computation  desired  followed 
by  vehicle  speed.  Use  a  1  for  user 
determined  routes  or  a  0  for  program 
determined  routes.  A  one  digit  entry  will 
designate  desired  speed  as  follows: 

1  for  9  mph 

2  for  12  mph 

3  for  15  mph 

4  for  18  mph 

Card  format  for  these  elements  is  12, 

IX,  12. 

X,  Y  coordinates  for  each  initial  red 
location.  One  card  is  needed  for  each 
red  unit  with  format  F6.1,  IX,  F6.1. 

If  the  user  has  selected  to  enter  his 
own  routes  in  data  element  5>  the 
following  is  required  for  each  route. 

a.  The  number  of  nodes  in  the  route 
(from  01. to  10).  This  card's  format 
is  12. 
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b.  X,  Y  coordinates  of  each  node 
along  the  route  in  order  from  v/est 
to  east.  Cne  card  for  each  node  in 
format  76.1,  IX,  76.1. 

3.  X,  Y  coordinates  for  blue  location,  force 

level  for  each  unit,  primary  search  directio 
for  that  unit  and  desired  search  width. 

Force  level  must  be  between  1.0  and  3«0> 
search  direction  between  135°  and  225°, 
and  search  width  from  30°  to  120°.  Cne 
card  for  each  blue  unit  with  format  76. 1, 

IX,  76. 1,  IX,  73.1,  IX,  13,  IX,  13. 

9.  Specify  if  blue  is  to  move  to  alternate 
defensive  positions,  breakpoint  distance, 
and  number  of  10  second  time  intervals 
allowed  for  the  move.  A  1  for  no  move 

or  0  for  option  to  move.  Format  for  this 
card  is  II,  IX,  76.1,  IX,  12. 

10. *  Alternate  blue  position  X  and  Y  coordinates 

if  move  is  desired  (data  element  9).  One 
card  for  each  location  with  format  76.1, 

IX,  F6.1 

11.  Fed  weapon  system  hit  and  kill  probab¬ 
ilities.  The  probabilities  are  entered 
for  6  range  bands,  with  one  card  for  each 
range  band.  Four  probabilities  for  each 

S"t 

range  are  probability  of  a  1  round  hit, 
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probability  of  a  hit  given  a  previous  round 
hit,  probability  of  a  hit  given  a  previous 


round  miss,  and  probability  of  a  kill 
given  a  hit.  The  6  range  bands  are  from 
0  to  500,  500  to  1000,  1000  to  1500, 

1500  to  2000,  2000  to  2500  and  2500  to  30C0. 
Each  of  these  six  cards  has  the  format 
4(F4.2,  IX). 

12.  Same  as  item  11,  except  the  probabilities 

for  the  blue  weapon  systems  are  entered. 

Six  cards  with  format  4(F4»2,  IX). 

In  the  formats  listed  above,  a  IX  is  a  space,  an  lx  is  an 
integer  of  x  digits,  and  a  Fx.y  is  a  real  number  of  length 
x  with  y  digits  to  the  right  of  the  decimal.  For  example, 
an  12  could  be  a  25,  an  F6.1  could  be  3487. 4»  A  format 
containing  these  specifications  could  be  12,  IX,  F6.1. 

The  data  card  would  then  look  like  25  3487.4,  with  the  25 
beginning  in  card  column  one. 
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card  1  2  3  4  ^?lj$ent 

column  1234567890123456789012345678901234567890  number 


1  28943 

1 

03  03 

2 

0000.0  2500.0  0500.0  4000.0 

3 

3.0  3.0  3.0 

4 

1  2 

5 

2500.0  7300.0 

2400.0  5900.0 

6 

2450.0  5300.0 

02 


4900.0  7200.0 
7100.0  6700.0 
02 

5350.0  6000.0  7 

7000. C  6300.0 
01 

6900.0  5500.0 

5700.0  7800.0  3.0  190  120 

5800.0  6500.0  3.0  180  120  8 

5900.0  5900.0  3.0  180  120 

0  0500.0  4  _  9 

7000.0  7200.0 

7100.0  6500.0  10 

7100.0  5800.0 _ 

0.85  0.85  0.75  0.70 
0.80  0.80  0.75  0.70 

0.75  0.75  0.70  0.65  11 

0.60  0.65  0.60  0.55 
0.45  0.50  0.50  O.35 

0.20  0.20  0.20  0.20  _ 


Figure  3«  Sample  Data  Deck 


card 

column 


1  3  4 

1234567890123456789012345678901234567890 


data 

element 

number 


0.60  0.70  0.65  0.85 
0.85  0.90  0.85  0.90 
0.80  0.85  0.85  0.80 
0.75  0.80  0.75  0.70 
0.60  0.70  0.65  0.65 
0.40  0.45  0.40  0.50 


Figure  3«  Sample  Data  Deck  (cont) 
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V.  REQUIRED  JCL 

The  model  and  line-of-sight  data  are  residing  on 
permanent  disk  in  the  W.  R.  Church  Computer  Center.  In 
order  to  exercise  the  model,  the  required  job  control 
language  (JCL)  is  illustrated  in  Figure  4.  Due  to  the 
CPU  time  required  for  execution,  a  time  parameter  of 
TIME  =  2  must  be  used  on  the  job  card. 

//JOB  CARD  ,  TIME=2 
//  EXEC  FORTCLG 

//FORT.SYSIII  DD  UNIT=2314, V0L=SER=PAT002,DISP=SHR , 

//  DS  N=  S 1 36  0  • SMO LER . PGM , DC B= ( RECFM=F B , LREC L=  8 0 , 3 LKS I ZE= 3 200 ) 
A 

//G0.FT08F001  DD  UNIT=2314 > V0L=SER=PAT002, DISP=SHR , 

//  DSW=S1360.LOSDATA,DCB=(RSCFM=FB, LRECL=80,BLKSIZE=3200 ) 

A 

//GO.SYSIN  DD  •» 

USER  DATA  DECK 


/* 


Figure  4.  Required  JCL 


VI.  EXPECTED  OUTPUT 


Once  the  model  has  been  executed,  the  user  can  expect 
the  following  output: 

1  -  Program  listing  of  the  model. 

2  -  A  summary  of  the  initial  battle  conditions, 

including  starting  locations  and  options  selected. 

3  -  A  summary  of  battle  conditions  after  each  10 

second  time  step,  including  unit  locations,  force 
levels,  unit  status,  percentage  lost  and  targets 
on  each  unit's  target  list.  The  status  will 
show  one  of  the  following: 

0  -  Unit  alive,  not  firing 

1  -  Unit  alive  and  firing 

2  -  Unit  killed 

3  -  Unit  moving 

4  -  A  statement  of  the  reason  for  battle  termination. 

A  sample  of  the  initial  battle  summary  is  at  Figure  5  and 
the  time  step  summary  is  at  Figure  6. 
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INITIAL  BATTLE  INFORMATION 


UNIT 

X 

Y 

FORCE  LEVEL 

1 

1900.0 

7800.0 

6.0 

2 

1700.0 

6500.0 

6.0 

3 

1700.0 

4900.0 

6.0 

4 

5800.0 

7700.0 

3.0 

5 

7400.0 

6300.0 

3.0 

6 

6200.0 

5150.0 

3.0 

ATTRITION  IS  STOCHASTIC 


ROUTES  DETERMINED  BY  USER 

RED  VEHICLE  SPEED  IS  15.0 

BREAKPOINT  DISTANCE  IS  1500.0 

BLUE  WILL  MOVE  TO  ALTERNATE  POSITIONS 
ALTERNATE  POSITIONS  ARE: 


UNIT 

X 

Y 

4 

9300. 

0  7300.0 

5 

9800. 

0  6200.0 

6 

9200. 

0  5800.0 

RED 

KILL 

PROBABILITIES 

RANGE 

P 

PHH 

PHM 

PKH 

500 

0.95 

0.95 

0.95 

1  .00 

1000 

0.80 

0.83 

0.81 

0.95 

1500 

0.65 

0.70 

0.70 

0.90 

2000 

0.40 

0.50 

0.55 

0.75 

2500 

0.30 

0.40 

0.40 

0.75 

3000 

0.20 

O.30 

0.25 

0.60 

BLUE  KILL  PROBABILITIES 

RANGE 

P 

PHH 

PHM 

PKH 

500 

0.80 

0.80 

0.80 

1  .00 

1000 

0.80 

0.80 

0.80 

0.95 

1500 

0.70 

0.70 

0.70 

0.90 

2000 

0.65 

0.65 

0.65 

0.85 

2500 

0.45 

0.55 

0.45 

0.70 

3000 

0.45 

0.55 

0.45 

0.60 

Figure  5.  Initial  Battle  Summary 


4? 


TIME  = 

UNIT 

480  SECONDS 

X  Y 

FORCE  LEVEL 

STATUS 

LOST- PC T 

1 

4933.2 

7193.2 

5.0 

0 

0.160 

2 

4832.3 

5799.9 

4.7 

0 

0.214 

3 

4902.6 

4579.8 

6.0 

0 

0.0 

4 

5800.0 

7700.0 

0.0 

2 

1.000 

5 

7400.0 

6300.0 

3.0 

1 

0.0 

6 

6200.0 

5150.0 

3.0 

3 

0.0 

Figure  6.  Time  Step  Summary 


TARGETS 


2 
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APPENDIX  B 

Program  Listing 


1 

2 
9 
8 

5 

6 
7 
a 
a 
to 

11 

12 

19 
1« 
IS 
ia 

17 

18 
ia 

20 
21 
22 
29 
21 
29 
28 

27 

28 
29 

90 

91 

92 
99 
98 

95 

96 

97 

98 

99 
80 
81 
82 
89 
88 
8S 
88 
87 
89 
89 
SO 


COMMON  /GAP1/  1  PRO  1 R  (6)  ,  1SECM0(6)  ,MVT0IR(6)  ,  X  (6)  ,T  (6)  ,3P0  (8) 
COMMON  /6AP2/  TA  (2) ,  T 1  (2) .  Th  (21 .  TM  (2) .  TF 1  (2> ,  TF2  121 .  TF3  (2) , 
19(2,6)  . PHM  (2,6)  .PHM  (2.6)  .  PKH  (2.6)  .  TF  (2) 

COMMON  /GAPS/  NBU.NAU.FL  (6)  ,F0  (6)  .  N01  (3)  .  X1C  (9. 200) .  TIC  (9,200) , 
110 I R (9.200) , AV5P. ISPD 

1 ,  IUSTAT  (6)  .11  (6)  .LOST  (6,6)  , VISFRA, V1SFRB.S1ZETK, 

19IZETM.NT  (6) , NF  (61  .SRF.01SMAX. 

1NL0SC  (6. 61  .  VISFR  (6. 6)  .  RM1NTK.RMXTK.RH  I  NTH,  RMXTW,  OP.  TONFR.  TNKFR, 
1PTT  (9.9)  ,RF. POR  (6.6)  .RPOA  (6.6)  .LOR  (6. 6)  .NR  (6)  .OFL  (6)  .POL  (6) 

COMMON  /GRP8/  TPOL  (6)  ,0LDQ(6,6)  .0(6.6) 

COMMON  /GRPS/  LOT  (6, 6)  .ROT  (6. 6) 

COMMON  /HILLS/  XC  (100)  .  TC  (100)  .PERK  (100)  ,SX  (100)  .ST  (100)  .RHO  (100) 
COMMON  /HILLS/  SCALE  (1001  .THORHO  (100)  .THOSCL  (100)  .BASE 
COMMON  /HILLS/  NH1LLS 

COMMON  /COVER/  CXC  (ISO)  . CTC  (ISO) , CPERK  ( 1  SO) , CPXX  ( 1  SO) , CPTT  (ISO) 

COMMON  /COVER/  CPXT  (ISO) .NOVELS 

COHHON  /COUNTR/  KH. KHH. KV. KN. KGRS. KELL. HINT 

COMMON  /GRID/  LST  (10. 101  .NHL  (10. 10)  .LlSTH  (850)  .KHREP  (100)  .KTREP 
COMMON  /GRID/  LSTC  (10. 10)  .NC  (10. 10)  .LlSTC  (800)  .KCREP  (ISO) 

COMMON  /GRP8/  ALPHA (6) 

COMHON  /GRP7/  XR  (6)  .  TA  (6)  .  1H0VE  (6) 

C 

C  INITIALIZATION. 

C 

CALL  OVFLOU 
BL-0.0 
AL-Q.O 
MP-0 

PAI-22. 0/7.0 
ZL*.  00001 
CALL  SETUP 
C 

C  AERO  TERRAIN  DATA  FOR  LINE  OF  SIGHT 

C  CHECK  FOR  STOCHASTIC  OR  DETERMINISTIC  ATTRITION 

C 

HERO  (5,130)  1TR1T.I3 
190  FORMAT  (It.  IX.  IS) 

DO  192  1-1.6 

CALL  RANDOM (IS, TRAN. 1) 

ALPHA  (I)  -  (-2.-TRAN--21  ♦  (2.  -TRAN*.  3) 

192  CONTINUE 
C 

C  ACAD  IN  NUMBER  OF  BLUE  ANO  RED  UNITS 
C 

READ  (5. 200)  NBU.NRU 
200  FORMAT  (12.  IX.  121 
C 

C  INITIALIZE  MEAPQN  SIZES 
C 

SIZETK-2.S 
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SIZETW-2.5 


C 

C  HERO  IN  EFFECTIVE  WEAPON  RANGES 
C 

READ (S, 102)  RMlNTK,  RMXTK, RMINTW, RMXTW 

102  FORMAT 1F6. 1.1X.F6. 1.1X.F6. I.1X.F6. 1. 1X1 
C 

C  INITIALIZE  PM, RF, TOWFR. TNKFR  AND  NOO 
C 

PM-.352 
RF-.S 
TOWFR-. 03 
TNKFR*. 1 
NOD*  2 

00  101  1*1. NRU 
N01  <11-125 
101  CONTINUE 
K-NRU-l 
L-NRU-NBU 
00  111  1-1. L 
II  (11-0 
111  CONTINUE 
C 

C  READ  IN  FORCE  LEVELS  OF  EACH  RED  UNIT 
C 

RE  AO  (5. 1031  (FL(l).l-l.NRUl 

103  F0flHArc3cF3.1,iXn 
C 

C  CHECK  FOR  TYPE  OF  ROUTE  0ETERMIN1TI0N 
C 

READ (5. 1061  IRTE.1SP0 

106  FORMAT  <11,  IX.  Ill 

IF  IISPO.EO. II  AVSP-9. 0 
lFUSPO.EO.il  OST-HO. 232 
IF  11SP0.E0.2J  AVSP-1 2. 0 
IF  IISPO.EO. 21  0ST-53.6U3 
IF  IISPO.EO.  31  AVSP-15.0 
IF  IISPO.EO. 31  OST-67. 053 
IF  IISPO.EO. HI  AVSP-18.0 
IF  IISPO.EO. Ml  0ST-80.  463 
C 

C  RE AO  IN  INITIAL  RED  LOCATIONS 
C 

00  6  1*1. NRU 

RE  AO  <5. 1071  XIC  <1 , 1 1 .  TIC  (1 , 1 1 

107  FORMAT (F6. 1. IX.F6. 1) 

6  CONTINUE 

IFIIRTE.EQ.il  GO  TO  108 
00  2  I-l.NRU 
00  2  >2.125 
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101  TlCtl.JI-TlClI.J-ll  *0ST-  (J-l) 

102  XICd.JI»XICd.J-II  ♦OST*  (J-l) 

109  101fld,J)-0 

109  2  CONTINUE 

10S  GO  TO  109 

108  100  CALL  MUTE 

107  109  5UMR0-0.0 

108  00  3  I-i.NRU 

109  FOdl-FLd) 

110  SUMRO-SUMfiO»FOd) 

in  xm-xica.ii 

112  Y(1)-T!C(1,1) 

113  MVTQlRd) -IQlRd.  1) 

119  SPOd)-AVSP 

US  lUSTATd)-0 

116  lPR0Ifidl-101Ad.il 

117  lSECNDdl-120 

118  NFd)-l 

119  lld)»l 

120  3  CONTINUE 

121  C 

122  C  RE AO  IN  INITIAL  BLUE  LOCATIONS  _ 

123  C 

129  SUNBO-O.O 

125  00  9  l-K.L 

126  HERO  (S.  1 091  Xd)  ,Td)  .FLd)  .IPAOlfid)  .ISECNDd) 

127  109  FORMAT  (F6.1,JX.F6.I,IX,F3.1.1X.I3.1X.I3) 

128  FOd)-FL(l) 

129  SUNBO-SUMBO-FO  d) 

130  MVTOIR  (I)  -0 

131  SP0(ll-0.0 

132  1USTAT  d I «0 

133  lMOVEdl-0 

139  9  CONTINUE 

135  C 

136  C  CHECK  FOR  ALTERNATE  BLUE  POSITIONS  AND  READ  IN  IF  WANTED 

197  C 

138  READ (S. 900)  1ALT. BREAK. ITEM 

139  900  FORMAT  dl,lXtF6.1,IX,12) 

190  IF  dALT.EO.l)  GO  TO  901 

191  00  902  I-K.L 

192  READ  IS,  107)  XAIlI.TAd) 

193  902  CONTINUE 

199  901  0ELT-1Q. 

195  TAdl-20. 

196  Tldl-8. 

197  THd)-6. 

198  TNd)-10. 

199  TF 1 d) -1 . 

ISO  TF2  d)  -I . 
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151 

TF3U1-1. 

152 

TA  (21-20. 

159 

T1  (2)  -8. 

154 

TH  (21  -8. 

155 

TH 121-15. 

156 

TF1  (21-10. 

157 

TF2121-12. 

158 

TF3121-15. 

159 

C 

160 

C 

REDO  IN  HIT  ANO  KILL  PROBABILITIES 

161 

C 

162 

DO  S  1-1.2 

163 

DO  514  J-1.6 

166 

REAO  (5.  Sl5J  P(l.  J1  .PHHd.Jl  ,  PHM  (I .  J1  .PKHC 

165 

515 

FORMAT  (4  (F4. 2,  1X1 1 

166 

514 

CONTINUE 

167 

5 

CONTINUE 

168 

PTT  it,  11-1.0 

169 

PTT  (1.21 -0.8 

170 

PTT  (2. 21  -0.2 

171 

PTT  (1,31  -0.8 

172 

PTT  (2,31-0.15 

173 

PTT  (3.31 -0.05 

174 

00  31  1-1. NRU 

175 

00  31  J-K.L 

176 

NLQSC(l.J)  -0 

177 

NLOSC  (J.  11 -0 

178 

Od.Jl-l.O 

179 

QU.  11-1.0 

180 

WISFR1I.  Jl-0.0 

18! 

VISFR  (J.  11  -0.0 

182 

31 

CONTINUE 

183 

1C-1 

184 

C 

18S 

C 

PRINT  INITIAL  BATTLE  INFORMATION 

186 

C 

187 

WRITE  (6.599) 

188 

599 

FORMAT  Cl  MX.  "INITIAL  BATTLE  INFORMATION 

189 

WRITE  16.600) 

190 

600 

FORMAT (/IX. "UNIT ",7X. "X'.8X. *T*,4X. "FORCE 

191 

DO  601  1-1. L 

192 

WRITE  (6,602)  I  ,X  (I)  ,T  (II  ,FL  (I) 

193 

602 

FORMAT (IX, I3.3X.F7.1.2X.F7.1.7X.F3. 1) 

194 

601 

CONTINUE 

195 

IF  (ITRIT.EQ.  1)  GO  TO  603 

196 

WRITE (6,604) 

197 

604 

FORMAT (/IX. "ATTRITION  IS  STOCHASTIC VI 

198 

60  TO  60S 

199 

603 

WRITE  (6,606) 

200 

606 

FORMAT (/IX. "ATTRITION  IS  DETERMINISTIC’/) 
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201  60S  IF  (1RTE.EQ.0)  GO  TO  60? 

202  MAI TE  (6,608) 

203  608  FORMAT  I1X, 'ROUTES  DETERMINED  BY  USER'/) 

204  607  WRITE  (6,609)  AVSP 

205  609  FORMAT  (IX. 'AEO  VEHICLE  SPEED  IS  *.FM.  1/1 

206  WRITE  (6,610)  BREAK 

207  610  FORMAT  (IX. 'BREAKPOINT  DISTANCE  IS  '.F6.1/I 

208  IF  (1ALT.EQ.Q)  GO  TO  61 S 

209  WRITE  (6,620) 

210  620  FORMAT  (IX. 'BLUE  WILL  NOT  MOVE  TO  ALTERNATE  POSITIONS’/) 

211  GO  TO  625 

212  615  WRITE  (6.630) 

213  630  FORMAT  (IX. 'BLUE  WILL  HOVE  TO  ALTERNATE  POSITIONS VIX. 

214  1  'ALTERNATE  POSITIONS  AREi ’/IX,  'UNI T ' . 5X. *X ' . 6X.  'T M 

215  DO  635  1-K.L 

216  WRITE  (6,640)  1 . XA  (I)  .  YA  (II 

217  640  FORMAT  (1X,13,3X,F7.1,2X,F7.1) 

218  635  CONTINUE 

219  625  IRAN-500 

220  WRITE  (6,645) 

221  64S  FORMAT  (/4X. 'RED  KILL  PR0BA81LITI ES  VI X.  'RANGE \4X.  'P ', 

222  14X,  ’PHH',3X,  'PHH".3X,  'PKH'I 

223  DO  650  1-1,6 

224  WRITE  (6,  655)  1  RAN.P  (1 .  11  .  PHH  (1 .  I)  .  PHM  (1 .  I)  ,PKH  (I .  I) 

225  655  FORMAT  (2X. 14.4  (2X.F4. 2)) 

226  1RAN-1RAN*500 

227  650  CONTINUE 

228  IRAN-500 

229  WRITE  (6.660) 

230  660  FORMAT  1/4X. '8LUE  KILL  PROBABILITIES  VI X.  'RANGE  ', 4X.  'P\ 

231  14X. 'PhH’,3X. 'PHM'.3X, 'PKH') 

232  00  665  1-1.6 

233  WRITE  (6.655)  IRAN. P  (2, 1)  , PHH  (2. 1)  .  PHM  (2, 1)  , PKH  (2, 1) 

234  1RAN-1RAN-*500 

235  665  CONTINUE 

236  WRITE  (6.670) 

237  670  FORMAT  Cl  MOX, 'BATTLE  BEGINS'//) 

236  C 

239  C  UPDATE  LOCATION  OF  RED  UNITS. 

240  C 

241  0ISMAX-5000.0 

242  67  DO  9  I-l.NRU 

243  IF  UUSTAT  (I )  . EO. 2)  GOTO  9 

244  IF  (1USTAT  (I)  .EQ.O)  GOTO  76 

245  NF  (I)  -NF  (I)  ♦! 

246  IF  (NF  (I)  . LT.NOD)  GOTO  S 

247  NF  (1)  -1 

246  76  00  11  J  -  1.  NRU 

249  IF  U  .EQ.  II  GO  TO  11 

250  IF  (lUSTAT(J)  .EQ.  2)  GO  TO  II 
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251 

252 

253 
25V 

255 

256 

257 
2S6 
256 
260 
261 
262 
263 
26V 

265 

266 

267 

268 

269 

270 

271 

272 

273 
27V 

275 

276 

277 

278 

279 

280 
281 
282 
283 
28V 

285 

286 
287 
268 
269 

290 

291 
282 
293 
29V 

295 

296 

297 

298 

299 

300 


OIST  -  X  (11  -  X(J1 
IF  (OIST  .GT.  OISHflXl  GO  TO  9 
11  CONTINUE 

11  (1)  -  11  (1)  *  1 
K7-1I  (11 
X  (11  -XIC  (1  ,K7) 

T  (1» -TIC  Cl .  K7I 
MVTOIR (I I  -I  DIR  (1  ,K7) 

IPROIRdl  •  1 01ft  (I.K7I 
9  CONTINUE 
C 

C  LINE— 0F-S1GHT  CHECK  BETWEEN  UNITS  AND  TARGETS  SELECTION 
C 

DO  17  J-K.L 
NT  (JJ  -0 
17  CONTINUE 

DO  12  1»1 . NRU 
NT  (II  >0 

IF  (1USTAT  (1)  .ED. 2)  GOTO  12 
DO  16  J-K.L 

IF  (IUSTAT  (J)  .  EQ.  2.  OR.  1  LIST  AT  ( J)  .  EQ.  3)  GO  TO  16 
XXl-X(l) 

TT1-T  (I) 

133  CALL  ELEV (XXI >  TT1 , THAC1) 

XX2-X  (J) 

TT2-T  (J1 

CALL  ELEV  (XX2, TT2. TMACJI 

LAT0B-1 

L8T0A-1 

CALL  LOS (XXI. TT1. TRACI. 0.0, S1ZETK, XX2, TT2.TNACJ.0.0.S1ZETH. 
1LAT08,LBT0A, VISFRA, VISFRB) 

VISFR  (1 ,  Jl  -VISFRA 
VISFR  (J.U-V1SFR8 
IF  (V1SFRA.GT.ZL)  GOTO  18 
LOST  (I.  JJ  -0 
LOST  (J.  1)  -0 

NLOSC  (I.  J)  -NLOSC  (I,  J)  *1 
NLOSC  (J.  11 -NLOSC  (I,  J) 

GOTO  16 
16  LOST  (1 ,  J)  >1 
LOST  (J. II -1 
NLOSC  Cl.JJ-0 
NLOSC  U.11-0 

RANGE-SORT  ( (X  (I)  -X  (JJ  J  (T  (I)  -T  (JJ )  «21 

IF  (RANGE. LT.RMINTK. OR. RANGE. GT.RHXTKJ GOTO  20 

IFIOII.JJ.EO.I.OJ  GOTO  20 

IUSTAT  (I) -1 

NT  (11  -NT  (11  -M 

N-NT  (11 

LOT  (1  ,NJ  -J 
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301 

ROT  (l.H)  •RANGE 

302 

IF  (H.EQ.l)  GOTO  20 

303 

CALL  SORT  (l.H) 

304 

20 

IF  (RANGE. LT.RH1NTH. OR. RANGE. GT.RHXTH)  GOTO 

16 

30S 

IF  <Q(J,II  .EQ.1.0)  GOTO  16 

306 

1USTAT  1J1  -1 

307 

NT  (Jl  -NT  (JJ  ♦  ! 

308 

M-NT  (Jl 

309 

LOT  U.H)  -1 

310 

ROT  U.H)  -RANGE 

311 

IF  (H.EQ.l)  GOTO  16 

312 

CALL  SORT  U.H) 

313 

16 

CONTINUE 

314 

12 

CONTINUE 

315 

00  25  1-1.  NRU 

316 

IF  (IUSTAT  (1)  . EQ. 2)  GOTO  25 

317 

IF  (NT  (I)  . NE.O)  GOTO  25 

318 

IUSTAT  (11  >0 

319 

NF  (l)-0 

320 

25 

CONTINUE 

321 

00  79  J-K.L 

322 

IF  (IUSTAT  (J)  .  EQ.  2.  OR.  IUSTAT  U)  .EQ.3)  GO 

TO 

79 

323 

IF  (NT  U)  .EO.O)  IUSTAT  (J) -0 

324 

79 

CONTINUE 

325 

C 

326 

c 

UPDATE  OF  THE  ACCUHULATEO  DETECTION  PROBABILITIES 

327 

c 

326 

IAA-1 

329 

IBB- NRU 

330 

1CC-K 

331 

IQQ-L 

332 

FR-TOWFR 

333 

OP-PM 

334 

37 

00  14  I-IAA.IBB 

335 

IF  (IUSTAT  (1)  .EQ. 2. OR.  IUSTAT  (I)  .EQ.3)  GO 

TO 

14 

336 

00  19  J-ICC.IOQ 

337 

PROP-Q.O 

338 

IF  (IUSTAT  (J)  .EO. 2. OR.  IUSTAT  (J)  .EQ.3)  GO 

TO 

19 

339 

OLOQ  (1 .  Jl  -0  (1 ,  J) 

340 

IF  (LOST  (I .  J)  .EO.O)  GOTO  15 

341 

IF(NTtl).GT.O)  GOTO  22 

342 

PCTVIS-VISFRIJ.ll 

343 

CALL  LAHOA  (I , J. PCT VI S. OETRAT , PSU8K) 

344 

OV-EXP  (-  (FL  (1)  «OETRAT»OP«DELT*FL  (J) ) ) 

345 

IF  (NT  (J)  . GT.O)  GOTO  23 

346 

fl(I.J)  -Q(I.J)  »QV 

347 

GOTO  19 

349 

23 

OP-  (1.0-PSU8K)  «>«  (FR-DELT-FL  (I)  «FL  (J)  1 

349 

0(1 .  J)  -0  (I .  J)  h  (OV+OP-QV-OP) 

350 

GOTO  19 

55 


951 

22 

N5»NT  11) 

352 

00  2V  11-1.M5 

953 

Kl-LOT  (1.111 

954 

ANG1-ATAN2  IT  (K 1 1  - T  ( 1 1  .  X  (Kl)  -X  11)  1 

355 

ANG2-ATAN2  (7  IJ)  -T  11)  ,  X  1J)  -X  11)  ) 

356 

IF  (1ANG1-ANG2)  .GE.O.O)  GOTO  77 

357 

IF  1ANG2.LT. 0.0)  GOTO  32 

3SB 

ANG*2*PA 1 ♦ ANG1 -ANG2 

3S9 

GOTO  35 

360 

32 

ANG*2»PA1 *ANG2-ANGl 

361 

35 

IF  (ANG.GT.PA1)  ANG*2«PAI-ANG 

362 

GOTO  33 

363 

77 

ANG*ABS  (ANG2-ANG1 ) 

36V 

33 

AA>1 5. OmPAI / 180. 0 

365 

IF  (ANG.GT.AA)  GOTO  2V 

366 

PROP«PROP*PTT (1 1 .NS) 

367 

2V 

CONTINUE 

368 

IF  (PROP. EQ. 0.0)  GOTO  3V 

369 

IF  (NT  1 J)  .  GT . 0)  GOT 0  36 

370 

CALL  LAMOAU. J.PCTV1S.0ETRAT.PSUBK) 

371 

DETRAT-OETRAT-RF 

372 

QV-EXP  1-  (FL  (1)  >«PHOP«OETRAT*DELT>.FL  U) ) ) 

373 

0(1.J).Q(1.J)«0V 

37V 

GOTO  19 

375 

36 

0(1,  J)  >0.0 

376 

GOTO  19 

377 

3V 

IFllAA.EO.il  GOTO  19 

378 

Oil.  J)»1.0 

379 

GOTO  19 

380 

15 

IF  (NLOSC  (1.  J)  .  LE .  3)  GOTO  19 

361 

mi.ji'i.o 

382 

19 

CONTINUE 

383 

IV 

CONTINUE 

38V 

IF  (1AA.EQ.K)  GOTO  38 

385 

FR-TNKFR 

386 

IAA-K 

387 

1BB-L 

388 

ICC-1 

389 

IOO-NRU 

390 

OP-1.0 

391 

GOTO  37 

392 

C 

393 

C 

FIRE  ALLOCATION. 

39V 

C 

395 

38 

00  28  I«1.L 

396 

28 

NACD-0 

397 

00  26  1*1 .L 

998 

IF  (IUSTAT  (II  . EO.  2.  OR.  IUSTAT  (I)  .EQ.  3)  GO  TO  26 

399 

IF  (NT  (1)  .  EQ.  0)  GOTO  26 

VOO 

DO  27  J-1,3 
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401 

402 

403 

404 

405 

406 
40? 

408 

409 

410 

411 

412 

413 

414 

415 

416 

417 
416 

419 

420 

421 

422 

423 

424 

425 

426 

427 

428 

429 

430 

431 

432 

433 

434 

435 

436 

437 

438 

439 

440 

441 

442 

443 

444 

445 

446 

447 

448 

449 

450 


APOA  (I.  J) -0.0 
27  CONTINUE 

IF  (NT  (1)  .£0. 1)  GOTO  78 
IF  (NT  (11  .E0.2I  GOTO  29 
NOT*  3 

MMl-tOT  (1.1) 

MH2-(.0T  (1.2) 

HM3-LQT  (1.3) 

PROB-  0.0-0  (1 .  MM1)  J  *0  (1 ,  MM2)  «Q  (I.  MM3) 

APOA  (1 . 1 )  -APOA  (1.1)  »PTT  (1 . ))  xPROB 
PROB-Q  (1  .KM1)  *  (1 .0-0  U  ,  MM2)  1  *Q  (1 .  MM3) 

APOA  (1.2)  -APOA  (1.2)  -PTT  (1.1)  xPROB 
PROB-Q  (1 ,  HH1 1  mO  (1 .  MM2)  h  (1 . 0-0  (1 .  MM3)  ) 

APOA  (1.3)  -APOA  (1.3)  -PTT  (1.1)  xPROB 

PROS*  1 1 . 0-0 ( I , MM1 ] ) *  ( 1 . 0-0  ( I . MM2) ) *Q ( 1 . HM3I 

APOA  (1. 1)  -APOA  (1,1)  -PTT  0.2)  xPROB 

APOA  (1.2)  -APOA  (1.2)  *PTT  (2.2)  xPROB 

PROS*  (1 . 0-Q  ( 1 .  MM  1 ) )  *Q  O  .  MM2)  *  O  .  0-Q  ( 1 .  MM3) ) 

APOA  (1. 1)  -APOA  (1.1)  -PTT  (1,2)  xPROB 
APOA  (1.3)  *APOA  (1.3)  *PTT  (2.2)  -PROB 
PROB-Q  ( 1 .  HH1 )  X  ( 1 . 0-0  ( 1 . MM2)  )  *  (1 . 0-0  1 1 ,  MM3) ) 

APOA  (1.2)  -APOA  (1.2)  *PTT  (1.2)  xPROB 
APOA  (1.3)  -APOA  (1.3)  -PTT  (2.2)  -PROS 
PROB-  (1 . 0-Q  (1 .  MH1)  )  -  (1 . 0-0  (1 ,  HH2) )  x  (1 . 0-0  (I .  MM3) ) 
APOA  (1 . 1)  -APOA  (1.1)  -PTT  0.3)  -PRQB 
APOA  (1.2)  -APOA  (1.2)  -PTT  (2. 3)  xPROB 
APOA  (1.3)  -APOA  (1.3)  *PTT  (3.3)  -PROB 
30  DO  44  J- 1 .NOT 
KK-LOT  (1.  J) 

NA  (KKJ-NA  IKK)  ♦  1 
IN-NA(KK) 

LOA  (KK,  IN)  -1 

POA  (KK,  1 N)  -  APOA  (I ,  J) 

44  CONTINUE 
GOTO  26 
29  NOT-2 

MM1-L0T  (1,1) 

HM2-L0T  (1 ,2) 

PROB- O . 0-0 ( 1 . MM 1 ) ) *Q O , MM2) 

APOA  (1.1)  -APOA  (1,1)  -PTT  (1,11  xPROB 
PROB-Q  (I. HMD  x  (1 .0-0  (I. MM2)  I 
APOA  (1.2)  -APOA  (1.2)  -PTT  (1.  II  -PROB 
PROS-  (1. 0-0(1. MM1)  )*  (1. 0-0  (I. MM2)  I 
APOA  (1.1)  -APOA  (1, 1)  +PTT  0,2)  -PROB 
APOA  (1.2)  -APOA  (1.2)  -PTT  (2.2)  -PRQ8 
GOTO  30 
78  NOT- I 

HM1-L0T  (1.1) 

PROB- 1 . 0-Q  (1 .  MM1) 

APOA  (1. 1)  -APOA  (1. 1)  -PTT  (1. 1)  xPROB 
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451 

GOTO  30 

452 

26 

CONTINUE 

453 

C 

454 

C 

ATTRITION. 

455 

C 

456 

SUMfl-0.0 

457 

SUM8-0.0 

458 

00  40  1-1. L 

459 

IF  (IUSTAT  ll)  .  EQ.2.0R.  1USTRT  111  .EQ.3)  GO  TO  40 

460 

H6-NAH1 

461 

SUM-0. 0 

462 

IF  (M6.  EQ.  0)  GOTO  47 

463 

00  41  J-1.M6 

464 

M?«L0A(1.  J» 

465 

1FIH7.LT.K1  GOTO  42 

466 

ITTPE-2 

467 

GOTO  43 

468 

42 

1TTPE-1 

469 

43 

RANGE-SORT  ( (X  111  -X  <H7)  1  -2-  (T  (1)  -T  IH7)  1  —2) 

470 

IF  (1TR1T.EQ.1)  GO  TO  1 31 

471 

CALL  STOCHIITTPE. RANGE. Ajll 

472 

GO  TO  5000 

473 

131 

CALL  ETKIITTPE. RANGE. T) 

474 

AJI-l.O/T 

475 

5000 

SUH-SUH-AJ1-FL  (M71  -PDA  (1.  J)  -DELT 

476 

41 

CONTINUE 

477 

47 

OFL  (11  -FL  (1) 

478 

FL  Ill  »FL  111  -SUM 

479 

IF  (FL  Ill  . GT.ZL1  GOTO  46 

480 

FL  (11-0.0 

481 

IUSTAT  (11  -2 

482 

46 

IF(l.LT.K)  GOTO  60 

483 

SUMB-SUHB-FL  (11 

484 

TPOL  (11  •  (FO  (I)  -FL  (I)  1  /FO  111 

485 

GO  TO  40 

486 

60 

SUMR-SUMR-FL (II 

487 

TPOL  (11-  (FO  (11 -FL  (11  1  /FO  dl 

488 

40 

CONTINUE 

489 

C 

490 

C 

PRINT  AND  CHECK  FOR  BATTLE  DETERMINATION. 

491 

C 

492 

IT IME- 1C- 10 

493 

00  57  I-K.L 

494 

IF  (IUSTAT  (II  .EQ. 21  GO  TO  57 

495 

00  56  J-l.NRU 

496 

IF  (IUSTAT  (J)  .EO. 21  GO  TO  58 

497 

CHECK-X  (1 1  -X  (J1 

498 

AVO-SORT  (  (X  (I)  -X  (Jl  1  mm2*  (T  (11  -T  (J)  1  mm2) 

499 

IF (AVO.LT. BREAK. OR. CHECK.LT. SO. 1  DO  TO  250 

500 

58 

CONTINUE 

58 


501  57  CONTINUE 

502  CO  TO  99 

SOS  c 

504  C  COMPLETE  BLUE  HOVE 

505  C 

506  250  00  251  1-K.L 

507  IF  (1  ALT.EQ.  1  .OR.  I  MOVE  (I)  .EQ.  ITEM)  CO  TO  6000 

SOB  IF  (IUSTAT  (11  .£0.0)  lUSTATdI-3 

509  lHQVEOl-lMOVEO)*! 

510  IF  (1H0VE  (I)  .LT.  ITEMI  CO  TO  2Sl 

511  Xdl-XA(I) 

512  Tdl-TAII) 

513  IF  UUSTAT  (I) . E0.3)  1USTATI11-0 

514  251  CONTINUE 

515  99  WRITE  16.112)  1 T 1  HE 

516  112  FORMAT C////1X. ’TIME  •  ’. 14. IX. 'SECONDS '//) 

517  WRITE  (6. 113) 

516  113  FORMAT  (IX,  'UNIT'.SX.  'X'.BX.  *Y’,5X,  'FORCE  LEVEL \2X.  'STATUS’. 

519  12X, 'L0ST-PCT'.2X. ’TARGETS'! 

520  00  59  I-l.L 

521  N6-NT  (11 

522  IF  (N6. NE. 01  GO  TO  4B 

523  WRITE  (6.2841  I . X  (11  .  T  (1)  ,FL  (II  .  1USTAT  (11  . TPOL  111 

524  264  FORMAT (3X. I 1.3X.F7. 1 . 2X. F7. 1 . 6X. F3. 1 . 9X. 1 1 . 6X. F5. 3) 

525  GO  TO  59 

526  4fl  WRITE  (6. 1141  1.  X  (11 .  T  111  .FL  UJ .  IUSTAT  ID  .TPOL  (J) , 

527  1  (LOT  (1.  Jl  .  J-1.N6) 

528  114  FORMAT  (3X.  1 1.3X.F7. 1.2X.F7. 1.6X.F3. 1 ,9X.  1 1 . 6X.F5. 3. 3X. 3  (1 1 ,  1X11 

529  59  CONTINUE 

530  C 

531  C  CHECK  FOR  BATTLE  DETERMINATION. 

532  C 

533  IOT-O 

534  00  53  I-l.NRU 

535  IF  (FL  (II  . EQ. 0. 01  GOTO  53 

536  JOT-1 

537  53  CONTINUE 

538  1F(I0T.EQ.U  GOTO  S4 

539  WRITE  (6.1 17) 

540  117  FORMAT  (IX. ’««•<  REO  FORCE  IS  ELIMINATED.  END  OF  BATTLE.  N 

541  GOTO  66 

542  54  I0T-0 

543  DO  55  1-K.L 

544  IF  (FL  (11  .EQ.O.OI  GOTO  55 

545  10T-1 

546  55  CONTINUE 

547  IFdOT.EQ.U  GOTO  65 

548  WRITE  (6. 118) 

549  118  FORMAT  (IX. 'h—  BLUE  FORCE  IS  ELIMINATED.  ENO  OF  BATTLE.’! 

550  GOTO  66 
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1  SUBROUTINE  LOStXft.TA.THACA.THlCR.SlZER.XB.YB.TNACB.THlCB.SIZQ. 

2  -LATQB.LBTOfi. V15FRA.  VJ  SFRB) 

s  c 

d  COMMON  /HILLS/  XC  (100)  .  TC  (100)  , PEAK  tlOO)  ,SX  (1001  ,ST (100)  .AHOUOQI 

S  COMMON  /HILLS/  SCALE  (100) .  TNORHO  (100)  .TW05CL  CQOl  .BASE 

0  COMMON  /HILLS/  NHlLLS 

7  COMMON  /COVER/  CXC  (ISO)  .CYC  (ISO)  .CPEAK  (150)  .CPXX  (ISO)  .CPTT  (ISO) 

B  COMMON  /COVER/  CPXT  (ISO) .NCVELS 

B  COMMON  /CflUNTA/KH,  KHM.  KV.  KN.  KGRS,  KELL.  K1MT 


10 

COMMON  /GRID/  LST  110.10)  .NHL  (10.10)  .LlSTH(dSO)  .KHREPdOOl  .KTREP 

11 

COMMON  /CR10/  LSTC  (10. 10)  .NC  (10.  >0)  .LlSTC  (400)  .KCREP  (ISO) 

12 

DIMENSION  1GX  (100)  ,  )GY  (100) .  IEL  (100)  .C31  (100)  .C32UOOI 

IS 

ORTA  NOR) O/iO/.GSIZE/ 1000./ 

Id 

C  SUBROUTINE  TO  COMPUTE  FRACTION  VISIBLE  FOR  OBSERVER  TARGET  PAIR3 

IS 

V1SFRA-I. 

IB 

V1SFR8-I. 

17 

XBA-XB-XA 

IB 

YBA-YB-TA 

19 

IF  ( (XBA.EO.O.)  .AND.  (YBA.EO.O.) )  RETURN 

20 

IFtSIZEA'THICA.LE.O.I  GO  TO  510 

21 

% 

IF  (S1ZEB*TH1CB.LE.0.)  go  to  5)0 

22 

IF (TM1CA.LT. 0.)  VI$FRA<1.0*TNICA/SIZEA 

29 

IF  ITN1CB.LT.0.)  V1SFRB-1 .0-TMICB/SIZEB 

29 

ZA-TMACA  ♦  TMICA  *  S1ZEA 

25 

ZB-TMACB  *  TMICB  «  SIZES 

20 

KTACP-KTRCP*) 

27 

ZBA-ZB-ZA 

2B 

XBASO<XBA<<XBA 

29 

YBASO<YBA<<TBA 

SO 

XTBA<XBA<<  TBA 

31 

TN0XBA<2. «XBA 

92 

THOTBA-2. «YBA 

S3 

C  COMPUTE  GRID  SQUARES  CROSSED  BY  A  TO  B  LINE 

3d 

NGRSO-O 

35 

IF(XBA)  110.9S.100 

SO 

95 

XBA<0. 1 

97 

100 

ISGX—1 

SO 

X1NC-0S1ZE/XBA 

99 

GO  TO  120 

dO 

110 

1SGX<1 

111 

XlNC— GSIZE/XBA 

«2 

120 

IF  (TBA)  IdO.  125.130 

d3 

I2S 

YBA-O. 1 

(Id 

130 

1SGT—1 

d5 

Y1NC-GS1ZE/TBA 

dO 

GO  TO  ISO 

07 

IdO 

1S6Y-1 

dB 

TINC— G31ZE/YBA 

09 

ISO 

1X<I*1FIX (XB/GSIZE1 

SO 

IF  (IX.GT.NGRIO)  1X-NGR10 

61 


SI 

IY»1*IF1X (YB/GS1ZE) 

52 

IF  (IY.GT.NGAIO)  1T-NGRID 

S3 

XNEXT-GS1ZE*  (FLOAT  (1X1  ♦0.5«  (1SGX-1.I) 

s« 

TNEXT-GS1ZE«  (FLOAT  (1Y)  ♦0.5*  (1SGY-1.1 1 

55 

XSTEP- (XB-XNEXT) /XBA 

ss 

YSTEP-  (YB-YNEXT)  /YBA 

57 

160 

NGRSO-NGRSQ*! 

58 

IGX  (NGRSO)-IX 

59 

1GT  (NGRS01-1Y 

80 

IFUXSTEP. GT.  1.1  .AND.  (YSTEP. GT.1.11  GO  TO  200 

6) 

IF  (XSTEP- YSTEP)  170.180,190 

82 

170 

1X-1X  +  ISGX 

83 

XSTEP»XSTEP*X1NC 

69 

GO  TO  160 

65 

160 

IX-IX+1SGX 

66 

XSTEP-XSTEP*X1NC 

67 

190 

IY«1Y*1SGT 

68 

YSTEP-YSTEP*YlNC 

69 

GO  TO  160 

70 

200 

KGRS-KGAS*NGRSQ 

71 

C  GRID  SQUARE  LISTJJ0N  COMPLETE  IN  IGX.  1GT  WITH  NGR3Q  ENTRIES 

72 

C 

73 

C  NON 

F1N0  WHICH  COVER  ELLIPSES  TOUCH  THE  A  TO  B  LINE. 

79 

C  CHECK  ELEVATIONS  AT  SI  ANO  S2  FOR  EACH  SUCH  ELLIPSE 

75 

NELS-0 

76 

CHTHAX-O. 

77 

IF  (NCVELS.EO-O)  GOTO  270 

78 

00  260  K-l.NGRSQ 

79 

1X-1GX  (K) 

80 

1Y-1GT  (K) 

81 

N-NC(IX.IY) 

82 

IF  (N.EO.O)  GO  TO  260 

83 

LS-LSTC(IX.IY) 

89 

LEND«LS*N-1 

as 

00  250  L-LS.LENQ 

86 

KELL-KELLM 

67 

IC-LISTC(L) 

88 

IF  (KCREPdCI  .EQ.KTREP)  GO  TO  250 

89 

KCREP(IC)  -KTREP 

90 

RX-XA-CXC(IC) 

91 

RT-YA-CTC(IC) 

92 

PXX-CPXX(IC) 

93 

PTT«CPTY (ICI 

99 

PXT-CPXY(IC) 

95 

AA>PXXxX6AS0*PrTHYBASa*PXYHXYBA 

96 

BB-PXX«THOXBAxRX«PYYi.THOYBA>«RT*PXY«  (RX«TBA«RYwXBftJ 

97 

CC»PXX«RX*<RX*PYY»«RTkRY«PXT«RX««RT-  1 . 0 

98 

ARG>B8«B6*9.Q«AA><CC 

99 

IF  (ARG.LE.O.)  GO  TO  2S0 

100 

SO-SORT  (ARG) 
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lot 

102 

109 

104 

105 
100 

107 

108 

109 

110 
111 
112 
US 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 


31—  (B8+SQ)  /  (2.0-AA) 

S2-  ISO-88)  /  (2.  OxAA) 

IFCSl.GE.l.)  CO  TO  250 
IF  (32. IE. 0.)  CO  TC  250 
IF (SI. IE. 0.1  GO  TO  S10 
IF (32. GE. 1.1  GO  TO  510 
C  NOW  CHECK  LOS  AT  SI  AND  S2 
K1NT-K1NTM 
CPK-CPEAKtlC) 

X5-XA-S2-XBA 

TS-TA«S2*TBA 

coll  euv(xs.rs.ms) 

HTS-MTS-CPK 

2S-ZA«S2«ZBA 

IF  (LAT08.EQ.01  GO  TO  210 

CALL  KOVER (ZA. TMACB, S I ZEB, ZB. S2. HTS. ZS. VI SFRBl 
IF  (VISFRB.LE.O.)  GO  TO  510 
210  IF  1LBT0A.EQ.0I  GO  TO  220 
3-1.-S2 

CALL  MOVER  (ZB,  THACA.  SI  ZEA.  ZA.  S.HTS.  ZS.  V1SFRA1 
IF  (VISFRA.LE.0.1  GO  TO  SlO 
220  XS-XA-Sl-XBA 
YS-TA-Sl-YBA 
CALL  ELEV  (XS.TS.HTS1 
HTS-HTS*CPK 
ZS-ZA*Sl«ZBA 

1F(LAT08.E0.01  GO  TO  230 

CALL  KOVER (ZA. TMACB. S1ZEB. ZB. SI . HTS. ZS. V1SFRB1 
IF (VISFRB.LE.O.)  GO  TQ  SlO 
230  lF(LBTOA.EQ.O)  GO  TO  240 
S-1.0-S1 

CALL  KOVER  (ZB. THACA.  SIZEA.  ZA.S. HTS.  ZS.VISFRA1 
IF  (VISFRA.LE.0.1  GO  TO  SlO 
240  NELS-NELSU 
I  EL  (NELS1  -IC 
CSl  (NELS)  -31 
CS2  (NELS)  -S2 

IFtCPK.GT.CHTMAX)  CHTMAX-CPK 
250  CONTINUE 
260  CONTINUE 
C  ACL  ECLIPSES  CHECKED 
C 

C  NON  START  ON  THE  HILLS 
270  00  600  K-l.NGRSO 
1X-IGX  (K) 

IT-IGT  (K) 

IF  (NHL  (IX.  IT)  .  EQ.Q)  GO  TO  600 
LS-LST  (IX.  IT) 

LENO-LS+NHL  (IX.  IT)  -1 
00  500  L-LS.LENO 
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151 

152 

153 
153 
155 
ISA 
157 

iss 

153 

ISO 

161 

162 

163 

IM 

165 

166 

167 

168 

169 

170 

171 

172 

173 
173 

175 

176 

177 

178 

179 
160 
161 
182 
163 
193 

185 

186 
187 
168 
168 

190 

191 

192 

193 
193 

195 

196 

197 
199 

199 

200 


I-LISTH(L) 

IF  (KHREP  (II  .EO.KTREPI  CO  TO  500 
KHREPdl-KTREP 

C  PROCESSING  FOR  HILL  1  STARTS  HERE 
KH-KH-1 

C  COMPUTE  M  -TOP  OF  THIS  HILL  ALONG  0-T  LINE 
CX-XBA/SX(I1 
CT-Y8A/ST  <11 
OX-  (XA-XC  (1)  1  /SX  (II 
OT-  (TA-TC  (II I  /ST  (I) 

FO-THOSCL  (II  «  <CX-OX«CT-OT-RHO  (1)  «  (CX*OY-CY«QX)  1 
GO-SCALE (11 « (CXmCX*CYkCY*TUQRHO (ii  *CX«CT1 
IF  (GO  .EO.O.l  GO  TO  500 
H— FQ/  (2.-G01 

IF  mes  (HI  .GT.S.l  GO  TO  500 
FSQ-FO-FO 

EQ-SCALE  111  -  (OXmOX*DThDT«TUORHQ(1)mDXhOT1 

PONER-EQ-FSQ/  (3.  -G01 

IF (POWER  .LT.  -3.1  GO  TO  500 

HHU-PEAK  (II  -EXP  (POUERl 

KHH-KHHK 

IF  (HHU.LE.BASE)  GO  TO  500 
ZH-ZA-HxZBA 

IFKN.LT. 0.1  .OR.  (N.GT.  1.11  GO  TO  300 
IF  (HHN  .GE.  ZHI  GO  TO  SlQ 
CVHTM-O. 

IF  (NELS.EQ.O)  GO  TO  300 
00  280  H-l.NELS 

IF  ( (CS1  (HI  .GE.H1  .OR.  (CS2  (HI  .LE.U)  1  GO  TO  2B0 
IC-IEL  (HI 

IF  (CVMTH.LT. CPEAK  (IC1  1  CVHTH-CPEAK  (IC1 
280  CONTINUE 

IF  ( (HHH-CVHTH)  . GE .  ZHI  GO  TO  510 
300  IF (HHH-CHTMBX.LT. AHIN1  IZA-SI ZEA, ZB-51 ZEBl  1  GO  TO  500 
C  IF  HE  GET  TO  HERE  THEN  NEED  TO  FIND  LOWEST  SIGHT  LINE  OVER  HILL 
C  NEWTON  ITERATION  A  TO  B  GIVING  VISFRB 
IF  (LATOB.EQ.  01  GO  TO  400 
KV-KV-1 
V-H 

HHV-HHH 

NCT-0 

FV-FO-V 

THOGV-Z.-GQ-V 

330  FCNV-ZA*HHVm  <TH0GV«V*FV-1.I 
KN-KN* 1 

FACTOR-  ITHQGVmTHOGV-2.  *  (GO^THOGV-FO)  -FSOl 

DFCNV-HHV-V-FACTOR 

IF  (ABS  (0FCNV1  .LT.l.E-101  GO  TO  3S0 

V-V-FCNV/OFCNV 

IF  (ABS  (VI  .GT.S.l  GO  TO  400 


FV-FQ-V 

TUOGV-2. -GQ-V 

POWER  •  EQ-F  V«GQ«V>*V 

IF  (POWER  .LT.  -3.)  GO  TO  100 

HHV-PEAK(l)  -EXP  (POWER) 

DHHV-HHV- (FQ+TWOGV) 

ELV-ZA-DHHV-V 

IF  (ABS  (HhV-ELV)  .LT.  1.)  GO  TO  350 
NCT-NCTM 

IF  (NCT.LT.  10)  GO  TO  330 
350  IF  ( (V.LT.O.) .OR.  (V.GT. 1.) )  GO  TO  100 
CVHTV-O. 

IF  (NELS.EO.O)  GO  TO  390 
00  380  H-l.NELS 

IF  ( (CSl  (H)  .GE.  VI  .OR.  (CS2  (HI  .LE.  V) )  GO  TO  380 
IC-IEL  (Ml 

IF  (CVHTV.  LT.  CPEAK  ( IC)  1  CVHTV-CPEAK  (20 
380  CONTINUE 
390  HTV-HHV»CVHTV 
ZV»ZA*V*ZBA 

CALL  KOVER (ZA.TMACB. SIZES. ZB. V.HTV, ZV. VlSFRB) 
IF  (VlSFRB. LE.O.)  GO  TO  SlO 
C  NEWTON  ITERATION  B  TO  A  GIVING  VISFRA 
400  IF  (LBTOA.EQ.O)  GO  TO  500 
KV-KV-1 
V«N 

VH1-V-1. 

HHV-HHN 

NCT-0 

FV-FQ-V 

TW06V-2.«G0»V 

430  FCNV-ZB^HHVk  (  (FO-»TWOGV)  «VM1  -1 . I 
KN-KN-1 

FACTOR-  (TW0GV-TH0GV+2.- (GO*TWOGV«FO) «FSO) 

OFCNV-HHVmVHI -FACTOR 

IF  (ABS  (0FCNV)  .LT.l.E-10)  GO  TO  450 

V-V-FCNV/OFCNV 

IF  (ABS  (VI  .GT.5.IG0  TO  500 

VH1-V-1. 

FV-FQkV 

TWOGV-2. mGQhV 

POWER  »  E0»FV«G0*<V«V 

IF  (POWER  .LT.  -3.1  GO  TO  500 

HMV-PEAK  (II  -EXP  (POWER) 

DHHV-HHV- (FO+TWOGVI 
ELV-ZB*OHHV«VHI 

IF  (ABS  (HHV-ELV)  .  LT.  I .  I  GO  TO  1150 
NCT-NCT*! 

IF  (NCT.LT . 1 0)  GO  TO  430 
150  IF  ( (V.LT. 0. 1  .OR. (V.GT. 1 . 1 1  GO  TO  500 


CVHTV-O. 

IF  (NELS.EQ.O)  GO  TO  490 
00  480  M-l.NELS 

IF  ((CS1  (Ml  .GE.V)  .OR.  (CS2  IH>  .LE.V)  )  GO  TO  480 
101EUH) 

IF  (CVHTV.LT . CPEflK  (IC1  )  CVHTV-CPERK  (IC) 

CONTINUE 

HTV-HHV«CVHTV 

ZV«Zfl*V*Z8fl 

S— VH1 

CALL  KOVER  (ZB. THflCfl. SI ZEfl, ZA.S.HTV.ZV, V  ISFRfl) 

IF  (V1SFRA.LE. 0. I  GO  TO  510 

CONTINUE 

CONTINUE 

RETURN 

VISFRfl-O. 

V1SFRB-0. 

RETURN 

ENO 
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m 


1 

2  C 

3 

4 

5 

6 

7 

8 
9 


10 

11 

12 

19 

47 

14 

7 

15 

17 

IS 

17 

18 

50 

19 

20 

21 

22 

23 

37 

24 

25 

26 

27 

28 

27 

29 

30 

60 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

C 

45 

100 

46 

47 

48 

49 

50 

SUBROUTINE  SETUP 

COMMON  /HILLS/  XC  (100)  ,TC(100)  .PERK  1100)  ,  SX  (100)  .ST  (100)  .  RHO  (100) 
COMMON  /HILLS/  SCALE ( 100) , TNORHO ( 1 00) . TWOSCL (100) . BASE 
COMMON  /HILLS/  NH1LLS 

COMMON  /COVER/  CXC  (150)  .  CTC  ( 150)  .  CPEAK  (150)  .  CPXX  (1  SO)  .  CPTT  (1 50) 

COMMON  /COVER/  CPXT (150) .NCVELS 

COMMON  /CQUNTR/KH, KHU, KV, KN, KGRS, KELL . K1 NT 

COMMON  /GRID/  LST  (10.  10)  .  NHL  (1 0.  1 0)  .  L ISTH  (ISO)  .  KHREP  (100)  .  KTREP 
COMMON  /GR10/  LSTCdO.  10). NC  (10.10). LISTC  (400). KCREP  (150) 

READ  18. 7)  NHILLS 
READ  (8. 47)  BASE 
FORMAT  IF  10.4) 

FORMAT  (16) 

FORMAT  (5F8.3.F6.4) 

00  50  1*1 . NHILLS 

RE  AO  (8.17)  XC  (I)  ,TC  (1)  .PEAK  (1)  .  SX  (1)  .  ST  (I )  ,  RHO  (1 ) 

CONTINUE 
READ  (8,37)  LST 
READ  (8.371  NHL 
READ  (8.7)  NHTOT 

RE  AO  (8. 37)  (L I STH  ( I )  .  1  - 1 .  NHTOT) 

FORMAT  (161  Si 
RE  AO  (8.71  NCVELS 
IF  (NCVELS. EQ.O)  GO  TO  65 
00  60  1*1 . NCVELS 

RE  AO  (8.27)  CXC  (I)  .CTC  (I)  .CPEAK  (1)  .CPXX  (1)  .CPTT  (1)  .CPXT  (I) 

FORMAT  (3F10. 4,  3E 1 3. 7) 

KCREP(l)—  2147483600 

CONTINUE 

RE  AO  (8.371LSTC 

REAO  (8.371  NC 

READ  (8.7)  NCTOT 

REAO  (8. 37)  (LISTC(l)  .1-1. NCTOT) 

65  00  100  1-1, NHILLS 
SX(I)-SX(I)xl.62S 
ST  (1)  -ST  (I)  -1.625 
XC  (I)  -  (XC  (1)  -500.1  -100. 

TC(l)-(YC(I)-930. )«100. 

TNORHO  (1)  -2. -RHO  (1) 

SCALE  ID—  l./(2.x  (l.-RHO(I)  xx2) ) 

TMOSCL  (I)  -2. -SCALE  (1) 

KHREP  (11  — 2147483600 

ALL  VALUES  NOH  IN  METERS  ON  0  —  10,000  GRID 
CONTINUE 

KTREP— 2147483600 

KH-0 

KHM-0 

KV-0 

KN-0 
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51 

52 

53 

54 

55 


KGRS'O 

KELL-Q 

KINT-0 

RETURN 

END 
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SUBROUTINE  ROUTE 


1 

2 

3 

4 

s 

6 

7 

8 

9 

10 
11 
12 
13 
19 
IS 
19 

17 

18 

19 

20 
21 
22 
23 
29 

25 

26 

27 

28 

29 

30 

31 

32 

33 
39 

35 

36 

37 

38 

39 

90 

91 

92 

93 
99 

95 

96 

97 
96 
99 
SO 


COHNON  /GRP3/  NBU.NRU.FL  (61 . FO  16)  .N01  (31  . XI C  (3.2001  . TIC  (3,2001 . 

UD1R  (3.200)  .flVSP.ISPO 

l.IUSTRT  16)  .11  (6)  .LOST  (6.6)  , VI 5FRA, V1SFRB, SIZETK, 

1 S 1 ZETH,  NT  (6)  .  NF  (61  .  SRF .  0 1  SHRX . 

1NL0SC  (6.6)  .V1SFR  (6.6)  .RMINTK.RMXTK.RMINTH.RMXTH.OP.TOHFR.TNKFR. 

1PTT  (3. 3)  ,  RF. POfl  (6.6). RPOR  (6.6). LOR  (6.6).  NR  (6)  , OFL  (6)  . POL  (6) 
DIMENSION  XLOC  (3.20)  .TLOC  (3.20)  . N  (3) 

IF  (1SPD.EQ.9)  DST-B0.963 
IF (ISP0.E0.3)  0ST-67.0S3 
IF (ISP0.EQ.2)  OST-53. 643 
IF  (1SPQ.EQ.  1)  0ST-40.232 
00  300  1-1. NRU 
RERO(S.IS)  Nil) 

IS  FORMAT  (12) 

NL-N(I1«1 

00  200  IN-2.NL 

REAO  (5.201)  XLOCS. TLOCS 

201  FORMAT (F6.1. IX. F6.1) 

XLOC  (I .  IN)  -XLOCS 
TLOC  (1.  IN) -TLOCS 

200  CONTINUE 

XLOC  (1. 1)  -XIC  (1, 1) 

TLOC  II.  1) -TIC  (I.  II 
10IR  (1. 1)  -0 
NL-N(l) 

NUH-2 

00  305  J-l.NL 

XL -XLOC  (I ,  J«  1 )  -XLOC  (I ,  J) 

TL-TLOC  (1 ,  J-* I)  -  TLOC  <1 ,  J) 

D1ST-S0RT  (XL--2*YL"«2) 

T-RBS  (TL) 

Z-T/XL 
RNGL-RTAN  (Z) 

OEC-RNGL -57.2958 
IFU.EO.l)  GO  TO  320 
XLN-  (OST-EXTRR)  -COS  (ANGL) 

01ST- (01 ST^EXTRR) -OST 
TLN- (OST-EXTRR) -SIN  (ANGL) 

XIC  (I.NUM)  -XIC  (I.NUM-1)  *XLN*XLE 
IF  (TL.GT.O.I  GO  TO  325 
TLN— TLN 

325  TIC  (I.NUM)  -TIC  (1  .NUM-l)  ♦TLN«TLE 
IF  (TL.GT.O.I  GO  TO  340 
I01R  (I.NUM)  — IFIX(OEG) 

GO  TO  341 

340  101R  (I.NUM) -IFIX(OEG) 

341  NUM-NUM*! 


69 


51 

320 

XLN-0ST«C0S ifiNGU 

52 

TLN-OSTnSIN (ANGL) 

53 

IF  (YL.GT.O.)  GO  TO  310 

54 

TLN— YLN 

55 

310 

lF(OlST.LT.DST)  GO  TO  315 

56 

XI C  (I , HUM)  »X1C  (1  , NUH-1)  ’XLN 

57 

TICU.NUHI-Y1C  (I.NUM-l)  »YLN 

sa 

IF  (YL.GT.O.)  GO  TO  342 

59 

lOlR(l.NUH)— IFIX(OEG) 

60 

GO  TO  343 

61 

342 

I01R  (1 ,  NUM)  >  IF  I X  (0EG1 

62 

343 

D1ST-01ST-0ST 

63 

NUM-NUM*1 

64 

GO  TO  310 

65 

315 

EXTRA-01ST 

66 

XLE-EXTRA«COS (ANGL) 

67 

YlE-EXTRAxSIN  IHNGL) 

66 

IF  (YL.GT.O.)  GO  TO  305 

69 

YLE— YLE 

70 

305 

CONTINUE 

71 

300 

CONTINUE 

72 

RETURN 

73 

ENO 
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SUBROUT  1 NE  LAMDA  1 1 , J. PCT V 1 S. DETRRT . PK» 

C  / 

C  SUBROUTINE  TO  COMPUTE  DETECTION  RATE  /(DETRRT)  OF  TARGET  J 

C  BT  OBSERVER  1  GIVEN  THE  VISIBLE  FRACTION  (PCTV1S). 

C 

COMMON  /GRP1/  1PRDI R  (6)  ,  I SECWQ  (6)  ,  MVTD1R  (6)  ,  X  (6)  .  T  (6)  ,  SPO  (6) 
TCFACT-1.0 
'  ZEROL-Q. 00001 
PA I -22. 0/7.0 

7  0-  (ISECHD  (I)  -PA1/180.Q)  /2.0 

BBS-  (1.0/  (2.0«  (SIN(O)  -D-COS  (Dll )) 

IF  (ABS(BBB)  .LT.ZER0L1  BBB-0.0 
AAA-  (-BBB)  -COS  ID) 

IF  (ABS(AAA)  .LT.ZEROU  AAR-0.0 
0TANG-ATAN2  <  (T  (J)  -T  (I) ) .  (X  (J)  -XII))) 

PD-1PR0IR (I) -PA1/ 180.0 
IF  MPO-OTANG)  .GE.O.O)  GOTO  1 
IF  IPO. IT. 0.0)  GOTO  9 
ANGLE-2-PAI*0TANG-P0 
GOTO  10 

9  ANGLE-2-PAI »P0-07ANG 
10  IF (ANGLE. GT. PA1)  ANGLE-2-PA1 -ANGLE 
GOTO  2 

1  ANGLE-ABS  IPO-OTANGI 

2  IF  (ANGLE. GT.OI  GOTO  3 
DUP*PD*0 
OLOH-PO-O 

ANGLFT-OTANG* (15.0-PAI/180.) 

IF  (ANGLFT.GT.OUP)  ANGLFT-OUP 
ANGLRT-OTANG-  (15. -PAI /1 80. ) 

IF  (ANGLRT.LT.OLOH)  ANGLRT-DLOH 

PK-BBB- (SIN (ANGLFT) -SIN (ANGLRT) ) ♦RAR« (ANGLFT-ANGLRT) 

IF  (PK.LT.O.O)  GOTO  3 
IF  (PK.GT .1.0)  GOTO  5 
GOTO  8 

3  PK-0.0 
DETRAT-0.0 
GOTO  6 

5  PK-l.Q 

8  RANGE-SORT  ( (X  (J)  -X  (I) )  —  2*  (T  (Jl  -T  (ll )  ""21 
Rfi-0. 001 -RANGE /PCT V I S 
T0ANC-ATAN2  ( (T  (Jl  -T  (J)  I ,  (X  (II  -X  (Jl )  I 
AO-MVTDIR (J) -PAI/180.0 
HORVEL-ABS  (SPO  (J)  -SIN  (TOANG-AD) ) 

OENOM-1 .953+TCFACT-  (0.5978*2.188-  (RR— 21  -0. 5038-H0RVELI 
IF  (OENOM.LE. ZEROL)  OENOM-ZEROL 
OETRAT-O. 003* 1 . 08B/0EN0M 
OETRAT-OETRAT-PK 

6  RETURN 
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I 


1  SUBROUTINE  ELEV(X.T.TMAC) 

2  C 

9  COMMON  /HILLS/  XC  (1001  .  TC  (1001 .  PERK  (1001  ,SX  (1001 , ST  (1001 .  RHO  (100) 

<1  COMMON  /HILLS/  SCALE  (100)  ,  TWORHO  (100)  ,  TWQSCL  (1001  , BASE 

5  COMMON  /HILLS/  NH1LLS 

6  COMMON  /GRIO/  LST  (10. 10)  .NHL  CIO.  101  .LISTH  (450)  .KHREP  (1001  .KTREP 

7  COMMON  /GRIO/  LSTC  (10. 10)  .NC  (10. 10)  .L1STC  (400)  .KCREP  (ISO) 

B  OATA  NGR10/10/. GS1ZE/1000. / 

9  C  FUNCTION  TO  COMPUTE  TERRAIN  ELEVATION  FOR  GIVEN  X.  T  COORDINATES. 


10 

ZHAX-BASE 

11 

IXal+lFlX (X/GS1ZE) 

12 

IF (1X.GT.NGR10)  1X-NGR10 

13 

IT-1-IF1X (T/GS1ZE) 

14 

IF  UT.GT.NGR10)  1T-NGRI0 

IS 

IF  (NHL  (IX,  IT)  .EQ.O)  GO  TO  ISO 

16 

LS-LST  (IX.  IT) 

17 

LENO-LS«NHL  (IX.  IT)  - 1 

18 

00  100  L-LS.LENO 

19 

I  «LISTH  (LI 

20 

OX-  (X-XC  (1)  1  /SX  (11 

21 

OXSO-OX-QX 

22 

IF  (QXSO  .GE.9.)  GO  TO  100 

23 

or-  (T-rcdii/sr  in 

24 

QTSO-QT-OT 

25 

IF  (QTSO  .GE.  9.1  GO  TO  100 

26 

0XT-TM0RH0(1)-0XmQT 

27 

FACTOR-SCALE  (11  M  (0XS0-*QTS0-»QXT) 

26 

IF  (FACTOR.LT. -3.)  GO  TO  100 

29 

HT-PEAK  (1)  -EXP  (FACTOR) 

30 

IF(HT.LE.ZMAX)  GO  TO  100 

31 

ZMRX  -  HT 

32 

100 

CONTINUE 

33 

ISO 

TNAC-ZMAX 

34 

RETURN 

3S 

END 
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1  SUBROUTINE  STOCK  <1 ,  RANGE ,  A1 

2  C 

3  C  SUBROUTINE  TO  CONFUTE  STOCHASTIC  ATTRITION  COEFFICIENT 

4  C 

5  COMMON  /GRP6/  ALPHA (6) 

6  COMMON  /GRP3/  N8U.NRU.Fl  (61  .FO  (61  .NOl  IS)  .XIC  (3.2001  .TIC  (3,200)  . 

7  1101R  (3,200)  .AVSP.ISPO 

8  l.IUSTAT  (61 .11  (6)  .LOST  (6.6) .  VISFRA,  VISFRB.SIZETK. 

9  1SIZETM. NT  (6)  ,NF(61  .SRF.01SHAX. 

10  1NL0SC  (6. 6) .  VISFR  (6. 6)  .  RHINTK.  RHXTK. RMl NTH.  RHXTH,  OP.  TOHFR.  TNKFR, 

11  1PTT  (3,3)  .RF.POA  (6.6)  .AP0A16.6)  .L0A(6.6)  .NA(6)  .OFL  (6)  .POL  (6) 

12  IF  (1  .EQ.2)  GO  TO  5003 

13  A-ALPHA  (1)  m  ( (1  .Q-RANGE/RHXTH)  ««2) 

14  GO  TO  5004 

15  5003  A-ALPHA (1)h(U.0-RANGE/RHXTK)hh2) 

16  5004  RETURN 

17  ENO 


1 

2 

a 

x 

5 

6 
7 
6 

9 

10 
11 
12 
19 

14 

15 

16 
17 
16 

19 

20 
21 
22 


SUBROUTINE  ETK(I.RANCE.T) 

c 

C  SUBROUTINE  to  compute  expected  time  to  kill  r  target. 

c 

COMMON  /GRP2/  TA  (21  ,  T1  (2)  .  TH  (2)  .  TM  (2)  .TF1  (2)  ,TF2  (2)  .  TF9  (2)  , 
IP  (2.61 . PHM  (2.61 ,  PHM  (2.6)  ,PKH  (2.6)  .TF  12) 

IF  (1  .EQ.21  GOTO  S 
TF  (I)  “TF  1  (!) 

GOTO  6 

5  IF  (RANGE. GT.  1000.0)  GOTO  7 

TF  (l)-TFl  (l)-(TFl  (Dm  (1000. 0-RANGE) /lQOO.Ol 
GOTO  6 

7  IF  (RANGE. GT. 2000.0)  GOTO  6 

TF  (I)  -TF2  (I)  -  ( (TF2  (II  -TF1  (!) )  *  (2000.0-RANGE)  / 1000.0) 

GOTO  6 

6  TF  (1)  -TF3  (1)  -  ( (TF3  (I)  -TF2  (II )»  (3000. 0-RANGE)  /1Q00.Q) 

6  J»  (RANGE*2S0.0) /500.0 

lFU.CT.6i  J«6 

T»TA  (1)  *T1  (1)  -TM  (1)  ♦  ( (TM  (I)  *TF  (II )  /PKH  (I .  J) )  ♦  ( (TM  (1)  *TF  (1)1/ 
1PHM  (I .  J)  I  m  ( (1 . 0-PhH  (I ,  J)  )  /PKH  II .  J)  *PHH  (1 .  J)  -P  (1 ,  J)  ) 

RETURN 

END 
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i 


c 


r 


1 

2  SUBROUTINE  SORT (1, HI 

3  COMMON  /GRPS/  LOT  (6.6)  .ROT  (6.6) 

4  DO  19  Jal.H 

5  IF  (ROT  (l.H)  . GE.ROT  (1.  J)  1  GOTO  19 

8  21  R-ROT(I.J) 

7  NN-LOT(l.J) 

8  ROT  (l.J) -ROT  (l.H) 

9  LOT  (l.J) -LOT  (l.H) 

10  ROT  (l.H)  -R 

11  LOT  (l.H)  -NN 

12  19  CONTINUE 

13  RETURN 

14  END 


i  __ 
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SUBRQUT IKE  KOVER  IZO.TMflCT . SI ZET ,ZT ,S,HTS, ZS, VlSFRT) 


1 

2 

9 

« 

5 

8 

7 

8 

9 

10 
11 
12 

13 

14 

15 


2000 

20S0 


IF  (S.EQ.O.I  GO  TO  2000 
IF  (HTS.GE.ZS)  CO  TO  2050 
HEXT-ZO- (HTS-ZO) /S 
EVIST-fiMAXl  IMEXT, THACT) 

IF  (EVlST.CE.ZT)  GO  TO  2050 
IF1EV1ST.LE.2T-SIZET)  RETURN 
VIS-  (ZT-EVIST)  /SIZET 
IF  (V1S.LT. VlSFRT)  VISFRT-V1S 
RETURN 

IF  (NTS. LT. 201  RETURN 

VISFRT-O.O 

RETURN 

ENO 
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APPENDIX  C 


Definition  of  Variables  in  Computer  Program 


ALPHA(I) 

APOA(I , J ) 

AVSP 

BREAK 

DISK AX 

DIET 

DST 

FL(  I ) 
FC(I) 

IALT 

IC 

IDIR(I,J) 
II  ( I ) 
IKOVE 

IPRDIR(I) 

IRTE 


Initial  attrition-rate  coefficient  for 
stochastic  attrition  module. 

±  Vi 

The  average  proportion  of  the  j  attacker 
of  unit  i  allocated  to  fire  on  unit  i. 

Average  speed  of  moving  attacking  units. 

Breakpoint  distance  between  attackers  and 
defenders. 

Maximum  distance  allowed  between  attacking 
units  before  the  leading  unit  is  delayed. 

The  straight-line  distance  between  two 
movement  nodes  inputed  by  the  user. 

The  distance  to  be  moved  each  time  step 
by  attacking  units. 

Force  level  of  unit  i. 

Initial  force  level  of  unit  i. 

Denotes  whether  the  user  desires  alternate 
defensive  positions  or  not. 

Counts  number  of  time  units  a  defender 
has  been  moving. 

x.  u  ^  Vi 

Direction  of  j  interval  in  i  route. 

Interval  index  for  unit  i. 

Number  of  time  units  a  defender  is  allowed 
for  moving  to  an  alternate  position. 

Primary  direction  of  movement  for  unit  i. 

Denotes  whether  user  wants  to  input  routes 
or  not. 

Random  number  seed  used  for  stochastic  attrition. 


IS 

I3EC'.7D(  I ) 


'.Vidth  of  search  sector  for  unit  i 


ISPD 

ITEM 

ITIME 

ITRIT 

IUSTAT(I) 
ICA(I,  J) 
LOST (I , J ) 

LOT  ( I ,  J ) 
MVTDIR(I) 
N(I ) 

NA(  I ) 

NBU 
inr  ( i ) 

NLCSCd,  J) 

NO 1(1) 

ITPvU 
NT  ( I ) 
OFL(I) 

P(I,J) 

PHH (I , J ) 


=  Input  variable  to  denote  user's  desired 
speed  for  attackers  movement. 

=  Input  variable  denoting  number  of  time 
steps  allowed  for  defender's  move. 

=  Current  time,  in  seconds,  of  battle. 

=  Input  variable  denoting  whether  attrition 
will  be  stochastic  or  deterministic. 

=  Current  status  of  unit  i. 

=  The  number  of  the  attacker  of  unit  i. 

=  Denotes  whether  line-of-sight  exists 
between  unit  i  and  j  or  not. 

=  The  number  of  the  target  of  unit  i. 

=  Movement  direction  of  unit  i. 

=  Number  of  nodes  inputed  by  user  for 
route  i. 

=  Number  of  attackers  of  unit  i. 

=  Number  of  blue  units. 

=  Number  of  time  units  unit  i  is  allowed  to 
fire  at  the  same  location. 

=  Number  of  continous  time  steps  that 

line-of-sight  does  not  exist  between  unit 
i  and  unit  j. 

=  Number  of  intervals  in  the  ith  route. 

=  Number  of  Red  Units. 

=  Number  of  targets  of  unit  i. 

=  Force  level  of  unit  i  during  previous 
time  step. 

s  t 

=  Probability  of  1  round  hit  by  unit  i 
in  range  band  j. 

=  Probability  of  a  hit  following  a  hit 
by  unit  i  in  range  band  j. 
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h::(i,  j) 


i  ,J) 


PCA( I , J ) 


?CL(  I ) 


PIT (I ,  J  ) 


2:  :i:rriv 


RCT( I ,J ) 

sizzt:: 

SIZST’Y 

TA(K) 


TF1(K) 


TF2(K) 


TF3(K) 


Probability  of  a  hit  follov/ing  a  miss 
by  unit  i  in  range  band  j. 

Probability  of  a  kill  given  a  hit  by 
unit  i  in  range  band  j. 

The  proportion  of  time  a  moving  unit 
is  searching  for  targets. 

The  proportion  of  the  attacker  of 
unit  i  allocated  to  fire  on  unit  i. 

Percent  of  unit  i  lost  during  the 
current  time  step. 

Proportion  of  surviving  firencwer 
allocated  to  the  i^h  target  if  there 
are  j  targets  available. 

Current  minimum  distance  betv/een  attackers 
and  defenders. 

Minimum  effective  range  for  attacking 
weapon  system. 

Minimum  effective  range  for  defending 
weapon  system. 

Maximum  effective  range  for  attacking 
weapon  system. 

Maximum  effective  range  for  defending 
weapon  system. 

The  range  of  the  jtn  target  of  unit  i. 

Size  of  attacking  vehicle. 

Size  of  defending  vehicle. 

Time  to  acquire  a  target  for  k  :  weapon 
system  type  (k  =  1,2). 

Tine  of  flight  to  1000m  for  k^^1  weapon 
system  type  (k  =  1,2). 

Time  of  flight  to  2000m  for  k^*1  weapon 
system  type  (k  =  1,2). 

Time  of  flight  to  3000m  for  kth  weapon 
system  type  (k  =  1,2). 


TH(K) 

TI(K) 

TI'(K) 

rpT[T<T?'nj 

T017FR 

TPOL(I) 

VISFR(I, J) 

VISFRA 

VISFRE 

:-:(!)  ,y(i  ) 

X  A  ( I ) ,  Y  A  ( I ) 

XIC(I,J) 

YIC (I , J) 

XL,YL 

XLCC(I, J) 
YLOC(I,J) 


=  Time  to  fire  a  round  following  a  hit 
for  weapon  system  type  k  (k  =  1,2)* 

=  Time  to  fire  first  round  after  target 
has  been  acquired  for  weapon  system 
type  k  (k  =  1,2). 

=  Time  to  fire  a  round  following  a  miss 
for  weapon  system  type  k  (k  =  1,2) . 

=  Firing  rate  for  attacking  weapon  system. 

=  Firing  rate  for  defending  weapon  system. 

=  Total  percentage  of  lost  since  battle 
began  for  unit  i. 

=  The  fraction  of  unit  i  seen  by  unit  j. 

=  Fraction  of  unit  A  as  seen  by  unit  E. 

=  Fraction  of  unit  B  as  seen  by  unit  A. 

=  Coordinates  of  unit  i. 

=  Coordinates  of  alternate  position  for 

defender  i. 

j.  \ 

=  Coordinates  of  the  j  interval  endpoint 
of  the  route  for  unit  i. 

=  Distance  added  to  previous  interval 
endpoint  for  vehicle  to  move  DST  during 
a  time  step. 

■t*  Vi 

=  Coordinates  of  the  j  node  inputed  by 
the  user  for  the  route  of  unit  i. 


80 


BIELICGRAPHY 


Bonder,  S. ,  "The  Lanchester  Attrition-Rate  Coefficient," 
Operations  Research,  v.  15,  p.  221-252,  May  196?. 

Bonder,  S.,  "The  Mean  Lanchester  Attrition  Rate," 
Operations  Research,  v.  18,  p.  179-181,  June  1970. 

Clark,  G.  M. ,  The  Combat  Analysis  Model,  Ph.D.  Thesis, 
Ohio  State  University,  Columbus,  Ohio,  1969. 

Karr,  A.  F.,  Stochastic  Attrition  Models  of  Lanchester 
Tijoe,  P-  1031,  Institute  for  Defense  Analysis,  Arlington, 
Virginia,  June  1974. 

Smoler,  J.,  Operational  Lanchester-Type  Model  of 
Small  Unit  Land  Combat,  M.S.  Thesis,  Laval  Postgraduate 
School,  Monterey,  California,  1979. 

Taylor,  J.  G.,  Force-on-Force  Attrition  Modelling, 
p.  20-89,  Military  Applications  Section  of  Operations 
Research  Society  of  America,  January  1980. 

Taylor,  J.  G.,  "Solving  Lanchester-Type  equations 
for  'Modern  '.Varfare'  with  Variable  Coefficients," 
Operations  Research,  v.  22,  p.  756-770,  October  1974. 

'.’/all ace ,  ’V.  S.,  and  Hagewood,  E.  G.,  Simulation  of 
Tactical  Alternative  Responses  (STAR).  M.  S.  Thesis, 
Laval  Postgraduate  School,  Monterey,  California, 

December  1978. 


81 


INITIAL  DISTRIBUTION  LIST 


1.  Defense  Technical  Information  Center 
Cameron  Station 

Alexandria,  Virginia  22314 

2.  Library,  Code  0142 
Naval  Postgraduate  School 
Monterey,  California  93940 

3.  Department  Chairman,  Code  55 
Department  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  California  93940 

4.  Professor  James  G.  Taylor,  Code  55Tw 
Department  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  California  93940 

5.  Professor  Glenn  F.  Lindsay,  Code  55Ls 
Department  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  California  93940 

6.  Professor  Samual  H.  Parry,  Code  55Py 
Department  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  California  93940 

7.  Captain  Glenn  M.  Mills 
710  Manhasset  Road 
Charlotte,  North  Carolina  28209 

8.  Command  and  General  Staff  College 
Attn:  Education  Advisor 

Room  123,  Bell  Hall 

Ft.  Leavenworth,  Kansas  66027 

9.  Headquarters,  Department  of  the  Army 
MILPERCEN 

Attn:  DAPC-OPF-C 

200  Stovall  Street 
Alexandria,  Virginia  22332 


82 


